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ABSTRACT 

Using the latest stellar evolution models, theoretical stellar spectra, and a compi- 
lation of observed emission line strengths from Wolf-Rayet (WR) stars, we construct 
evolutionary synthesis models for young starbursts. We explicitly distinguish between 
the various WR subtypes (WN, WC, WO), whose relative frequency is a strong function 
of metallicity, and we treat O and Of stars separately. 

We calculate the numbers of O and WR stars produced during a starburst and 
provide detailed predictions of UV and optical emission line strengths for both the WR 
stellar lines and the major nebular hydrogen and helium emission lines, as a function of 
several input parameters related to the starburst episode. We also derive the theoretical 
frequency of WR-rich starbursts. 

Our models predict that nebular He ll A4686 emission from a low-metallicity star- 
burst should be associated with the presence of WC/WO stars and/or hot WN stars 
evolving to become WC/WO stars. In addition, WR stars contribute to broad com- 
ponents beneath the nebular Balmer lines; the broad WR component may constitute 
several percent of the total flux in the line. 

We review the various techniques used to derive the WR and O star content from 
integrated spectra, assess their accuracy, and propose two new formulae to estimate the 
WR/0 number ratio from UV or optical spectra. 

We also explore the implications of the formation of WR stars through mass trans- 
fer in close binary systems in instantaneous bursts. While the formation of WR stars 
through Roche lobe overflow prolongs the WR dominated phase, there are clear obser- 
vational signatures which allow the phases in which WR stars are formed predominantly 
through the single or the binary star channels to be distinguished. In particular at low 
metallicities, when massive close binaries contribute significantly to the formation of 
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WR stars, the binary-dominated phase is expected to occur at ages corresponding to 
relatively low H/9 equivalent widths. 

The observational features predicted by our models allow a detailed quantitative de- 
termination of the massive star population in a starburst region (particularly in so-called 
"WR galaxies" ) from its integrated spectrum and provide a means of deriving the burst 
properties (e.g., duration, age) and the parameters of the initial mass function of young 
starbursts. The model predictions should provide the most reliable determinations to 
date. They can also be used to test the current theories of massive star evolution and 
atmospheres and investigate the variation in stellar properties with metallicity. 

Subject headings: galaxies: starburst — galaxies: stellar content — H II regions — 
stars: Wolf-Rayet 
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1. Introduction 

Wolf-Rayet (WR) galaxies are a subset of 
galaxies in which broad emission lines from WR 
stars are detected in the integrated galaxy spec- 
tra (Conti 1991). (A recent review of the prop- 
erties of WR stars can be found in Maeder Sz 
Conti 1994.) Typically, lines of He (He IIA 
4686), C (C III/C IV A4650, C IVA5808), and 
N (N IIIA4640) are seen. The strongest features 
are usually those at A 4650-4690, whose blended 
emission is often referred to as the "WR bump" . 
The strength of the various stellar emission lines 
has been used to estimate the number of WR 
stars present in these galaxies. Typical values are 
in the range N{WR) ~ 100 -10^ (see e.g., Vacca 
k Conti 1992, hereafter referred to as VC92). 
Because WR stars are the short-lived descen- 
dants of the most massive stars (M ^ 35Mq), 
the detection of WR emission lines in a star- 
burst galaxy spectrum immediately places pow- 
erful constraints on several parameters charac- 
terizing the starburst episode: (1) the number of 
WR stars relative to the number of O stars must 
be large, and therefore the burst episode must 
have been short; (2) the initial mass function 
during the star formation episode must have ex- 
tended to large masses (at least beyond 35 Mq); 
and (3) the time elapsed since the burst ended 
must be less than a few Myrs. Thus, the pres- 
ence of large numbers of WR stars in a starburst 
galaxy can be used as a signpost indicating a re- 
cent burst of massive star formation. 

The first models which attempted to quan- 
tify these constraints on the WR populations in 
starbursts were introduced by Arnault, Kunth Sz 
Schild (1989). The comparison of their models 
with observations of 45 extragalactic H II re- 
gions and/or galaxies with detected WR emission 
lines led to several important conclusions which 
were largely confirmed by later studies: (1) The 
observed decrease of the WR bump/H/3 inten- 
sity ratio with metallicity (Z) and the observed 
variations in the WR and O star populations in 



the local environment are due to the same ef- 
fect, namely the strong influence of Z on stel- 
lar mass loss leads to a significant decrease in 
the WR population at low Z (Maeder, Lequeux 
& Azzopardi 1980); and (2) The derived WR/0 
number ratio indicates that star formation must 
have occurred over only a short period compared 
to the lifetime of massive stars. 

Synthesis models including more recent evo- 
lutionary tracks and predicting a large variety 
of observational quantities were constructed by 
Mas-Hesse &: Kunth (1991a), and later updated 
by Cervifio & Mas-Hesse (1994). The WR fea- 
tures included in their models are the 4650 WR 
bump and C IV A 5808. Mas-Hesse & Kunth 
(1991b, 1997) used these models to analyze the 
spectra of nine WR galaxies. 

The use of the WR bump to determine the 
WR population presents several difficulties. The 
feature is composed of several stellar emission 
lines, whose strengths can vary substantially de- 
pending on the relative numbers of the various 
types of WR stars present. Furthermore, the stel- 
lar lines are often blended with nearby nebular 
emission lines of Fe, He, or Ar. Kriiger et al. 
(1992) were the first to synthesize the separate 
He II and C components of the WR bump, in an 
attempt to circumvent the problem of nebular 
line contamination, and to determine the domi- 
nant WR subtype present. 

The effect of varying star formation rates, ini- 
tial mass function, age and metallicity on massive 
star populations was studied in depth by Meynet 
(1995), who first stressed the importance of WR 
stars of the WC sequence. Although the exten- 
sive models of Leitherer & Heckman (1995) in- 
clude a detailed treatment of WR and O stars, 
observational features due to WR stars are not 
predicted. The same also holds for the evolu- 
tionary models of Garcia- Vargas, Bressan & Diaz 
(1995). 

Discoveries (usually serendipitous) of WR fea- 
tures in extragalactic objects have resulted from 
studies covering a wide range of topics, from 
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the primordial He abundance determination (cf. 
Kunth & Sargent 1983, Kunth & Joubert 1985, 
Izotov, Thuan & Lipovetsky 1994, 1997a; Izo- 
tov & Thuan 1997), the nature of Seyfert galax- 
ies (Heckman et al. 1997), to starbursts in cool- 
ing flows (Allen 1995). A considerable number 
of new observations can be expected with the 
new generation of 8-lOm class telescopes. To 
date the total number of known WR galaxies 
is ~ 90 37 of which were included in the first 
catalogue compiled by Conti (1991). However, 
quantitative analyses of the WR and O star con- 
tent has been carried out for only few 20) 
objects for which high quality spectra are avail- 
able (see VC92, Kriiger et al. 1992, Mas-Hesse 
k Kunth 1991b, 1997; Meynet 1995; Schaerer 
1996c; Schaerer et al. 1997). 

In the present paper we construct new synthe- 
sis models (building upon the models of Schaerer 
1996c) tailored to the analysis of massive star 
populations in young starburst galaxies. Us- 
ing up-to-date input physics (stellar evolution 
tracks, stellar atmospheres, and a new compi- 
lation of empirical WR line fluxes), we provide 
detailed predictions of many stellar and nebu- 
lar features from the UV to the optical, most of 
which have not been included in previous models. 
By combining our new models with the observed 
strength of the WR emission features, as mea- 
sured with medium to high spectral resolution 
spectra of sufficient quality (typically S/N ;^ 20), 
one can place constraints on the burst parame- 
ters which are stricter than those derived from 
earlier models. With these latest models it is 
possible to constrain the IMF slope (cf. Schaerer 
1996c), study the WC and WN populations in 
starbursts (cf. Meynet 1995; Schaerer & Vacca 
1996; Schaerer et al. 1997) and to address the 
question of the nature of nebular He II emis- 
sion in extragalactic H II regions on quantita- 
tive grounds (Schaerer 1996c, Schaerer & Vacca 
1996). More generally the present models can 
be applied to determine accurate burst parame- 
ters (age, duration, IMF etc.), and to study mas- 



sive star populations in different environments 
and their influence on the surrounding interstel- 
lar medium. 

The paper is structured as follows: The model 
ingredients, the input parameters, and the syn- 
thesized quantities are described in §^. In §^ 
we present predictions of the O star and WR 
stars populations, including distinct WR sub- 
types, and discuss the expected frequency of WR- 
rich starbursts. §Q summarizes issues related to 
the evolution of the ionizing spectrum of burst 
models. The predicted evolution of nebular H 
and He lines is discussed in §^; the evolution of 
WR lines is discussed in §^. In ^ we explore the 
effect of massive close binary systems on the for- 
mation of WR stars in young bursts. Different 
methods to derive the WR and O star content 
from integrated UV or optical spectra are pre- 
sented and re-discussed in §^. In §P we discuss 
applications and tests of our models and briefly 
sketch future prospects for studies of massive star 
populations in starbursts. 

2. Evolutionary synthesis models 

In this Section we describe the adopted model 
ingredients for our evolutionary synthesis mod- 
els, the most important input parameters, and 
the synthesized quantities. 

2.1. Evolutionary tracks, stellar spectra, 
and nebular continuum 

As in the case of most synthesis models we as- 
sume that for most purposes stellar evolution can 
be sufficiently well described by evolution models 
for single stars. The single star assumption forms 
the basis of what we call our "standard model". 
However, we also present calculations which in- 
clude close massive binaries in an approximate 
fashion. This first order approach will allow us 
to investigate the uncertainties in the results of 
the standard model arising from the neglect of 
binaries. 
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2.1.1. Single star evolution 

We use the recent evolutionary tracks of the 
Geneva group, which incorporate five different 
metallicities covering the range from Z=0.001 
(1/20 Zq) to Z=0.04 (2 Zq) (see Meynet et 
al. 1994 and references therein). In particular, for 
massive stars we use the recent grids of Meynet 
et al. , which incorporate mass loss rates that are 
a factor of two higher than those adopted in pre- 
vious evolutionary models. These high-mass-loss 
models reproduce a large number of statistical 
observations of massive stars, including the WR 
to O star ratios in regions of constant formation 
(Maeder Sz Meynet 1994) and they should there- 
fore be most appropriate for the purposes of the 
present work. 

Previous evolutionary synthesis models which 
explicitly incorporated WR stars used different 
evolutionary tracks. For example, the models of 
Cervifio & Mas-Hesse (1994) employ the Geneva 
tracks with standard mass loss rates, which yield 
smaller WR populations than the preferred high- 
mass-loss models. Leitherer Sz Heckman (1995) 
adopted older tracks from Maeder (1991b), for 
which detailed comparisons with observations are 
shown in Maeder (1991a). Garcia- Vargas, Bres- 
san & Diaz (1995) use the Padua evolutionary 
tracks for which no detailed comparisons regard- 
ing individual WR stars and WR populations 
have been published. We feel that the high-mass- 
loss Geneva models used in this paper should pro- 
vide the most accurate predictions of the WR and 
O star populations in a starburst region. 

2.1.2. Evolution of close massive binaries 

It is fairly well established that at low metal- 
licity the primary mechanism for the formation 
of WR stars is through mass transfer during the 
evolution of massive binaries, a formation chan- 
nel which should produce primarily early-type 
WN stars (Maeder 1982, Maeder 1991a, Maeder 
& Meynet 1994). Using the same evolutionary 
tracks as adopted in this work, Maeder & Meynet 



(1994) have shown that for regions of constant 
star formation with metallicities between ~ 1/10 
Zq (SMC) and ~ 1.75 Zq (M31), both the ob- 
served ratio of the total number of WR stars to 
the number of O stars present (the WR/0 star 
ratio) and the subtype distribution of WR stars 
(WC/WR, WNL/WR, and WNE/WR) can be 
well explained if a few percent of the O stars 
become WR stars of mainly WNE subtype as a 
consequence of Roche lobe overflow (RLOF) . 

To explore the implications of WR stars formed 
through mass transfer in massive close binaries 
we proceed in the following simplified manner. 
We adopt the recent evolutionary calculations 
from de Loore & Vanbeveren (1994) at differ- 
ent metallicities assuming an initial mass ratio 
of q =0.6, and Case B mass transfer. The cal- 
culations of de Loore & Vanbeveren (1994) in- 
clude initial primary masses from 9 to 40 Mq 
and neglect the possibility of subsequent WR 
formation of the secondary. Primaries with ini- 
tial masses Mi ^ 40-50 Mq should in general 
avoid Roche lobe overflow (RLOF) (cf. Vanbev- 
eren 1995); even if they do indeed experience 
RLOF their evolution should be nearly indistin- 
guishable from that of single stars (Langer 1995). 
For initial masses Mini ^ 9 Mq the formation 
of WR stars is not expected. Since the calcula- 
tions of de Loore &: Vanbeveren (1994) are not 
available at exactly the same metallicities as the 
single star models, we used the tracks with metal- 
licities closest to those for the single star models. 
(We adopted their solar metallicity models for 
Z >0.02, their LMC models for Z =0.008, and 
their SMC models for Z <0.004). 

Our primary goal is to derive, as a function of 
the age of the burst, the total number of WR 
stars (and their subtype distribution) formed 
from mass transfer in binary systems. For this 
purpose it is necessary to determine for each 
initial mass, whether or not at the given burst 
age the primary has already experienced RLOF 
leading to WR formation, and if it has, what 
the resulting WR subtype is. This is done in a 
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straightforward manner from the lifetimes in the 
WR phases as given by de Loore & Vanbeveren 

(1994) 

One additional parameter / is required to 
derive the total number of WR stars formed 
through the binary channel, where / is the ini- 
tial fraction of stars that are primaries in close 
binary systems (described by the parameters 
given above) and will thus experience RLOF dur- 
ing their evolution. Because the range in ini- 
tial masses covered by the binary models cor- 
responds to zero-age main sequence (ZAMS) O 
stars, our definition of / turns out to be essen- 
tially identical to the one given by Vanbeveren 

(1995) . In contrast Maeder k Meynet (1994) use 
^ = WRcb/0, which represents the fraction of 
O stars which become WR stars through RLOF. 
Note that both definitions implicitly neglect the 
evolution of the secondary, which would in prin- 
ciple affect the number of both O and WR stars. 
In the case of a constant star formation rate ^ 
and / are related by 

* ^ / [t(WRcb)/t(Oeb)]' 

I^Z7t{WR,^)/t{0MM)dM 

where t(WRcb) and t(Ocb) are the lifetimes of 
the primary of initial mass M in the WR phase 
and in the preceding O star phase, respectively. 
$(M) is the initial mass function, and Mmm 
and Mjnax are the minimum and maximum ini- 
tial mass of primary stars leading to WR for- 



^We note that the H-burning hfetimes tn from de Loore & 
Vanbeveren (1994, Table 1) are larger (typically by ~ 20 
%) than those from the single star models used here. The 
O star lifetimes given in Vanbeveren (1995) agree, how- 
ever, with those of standard evolutionary models (Schaller 
et al. 1992). To obtain a more consistent time for the on- 
set of Case B Roche lobe overflow we therefore adopt the 
in lifetimes from Meynet et al. This leads to an earlier 
appearance (~ 1 Myr) of WR stars formed from mass 
transfer than would be obtained by using the lifetimes 
from de Loore & Vanbeveren. The duration of the binary 
WR phases are taken directly from de Loore & Vanbev- 
eren (Table 1). 



mation by RLOF. To first order the models of 
Vanbeveren (1995) give mass- independent values 
for t(WRcb)A(Ocb) on the order of 0.2 and 0.1 
for Z=0.02 and Z < 0.01 (see his Figs. 2a-4a), 
and hence ^' ^ (0.2 - 0.1) x /. 

For numerical comparisons with Vanbeveren 
(1995) the standard mass loss tracks are the most 
appropriate. As Fig. 12 of Maeder & Meynet 
(1994) shows, the value of ^' required to repro- 
duce the observed WR/0 number ratio ~ 
0.08 to 0.01 for Z < 0.02) is fully compatible with 
the values of / advocated by Vanbeveren (1995). 
Yet another definition of the fraction of massive 
close binaries is adopted in the detailed models 
of Vanbeveren, Van Bever & de Bonder (1997). 
From their Table 1 (solar metallicity models us- 
ing observational constraints) one can, however, 
easily derive the values resulting for ^. One finds 
^' ~ 0.06-0.08, as expected for standard evolu- 
tionary models (see Maeder & Meynet, Fig. 12). 
This indicates that the assumed binary frequency 
of Vanbeveren, Van Bever & De Bonder is in 
agreement with the work of Maeder & Meynet 
(1994). From the comparison between the the- 
oretical and observed WR/0 ratio in regions of 
constant star formation using the high mass loss 
tracks (Fig. 11 of Maeder & Meynet 1994) and 
the approximate (metallicity dependent) conver- 
sion of \I' to / (see above) one obtains / ~ 0.1 
to 0.25 for Z=0.002 to 0.02. For the calculations 
presented here we adopt / = 0.2. 

Knowing the number of WR stars formed by 
mass transfer, we now proceed to derive their 
impact on the most important observable quan- 
tities, i.e. on nebular recombination lines and on 
the emission in the broad WR features. To this 
end one only requires a knowledge of the ion- 
izing flux and the number of WR stars in the 
different WR phases (cf. § [2.2.11 ). We adopt an 
average Lyman continuum flux of QoiyVK) = 
10^^ photons s~^, independently of the WR sub- 
type. This value roughly corresponds to the av- 
erage contribution of (single) WR stars derived 
from burst models at the time (~ 5 Myr) where 



6 



binary stars are first expected to be formed (cf. 
Fig. |A6| ). As a comparison the value used by 
VC92 is Qo(WR) = 1.7 X 10^^^ The contribu- 
tion of the ionizing flux in particular affects the 
H recombination lines, which are generally used 
as a reference for relative line intensity measure- 
ments. The broad WR line emission is treated 
as for single stars (see § 2.2.1). Since no infor- 



mation is available on the WC subtype from the 
binary models they are assigned the WC4 sub- 
type corresponding to the dominant subtype at 
low Z where the role of binaries is of largest im- 
portance. 

Our procedure completely neglects the de- 
tailed evolution of the binary stars in the H-R 
diagram. Except for the Lyman continuum, we 
neglect the continuum flux contribution of pri- 
maries in close binary systems, which affects only 
the predicted equivalent widths. This assump- 
tion is reasonable as long as the UV and optical 
continuum fluxes are dominated by single stars or 
stars in wide binary systems. The secondaries are 
assumed to evolve as single stars. The expected 
increase of their Lyman continuum flux after 
mass transfer (due to possible "rejuvenation") 
can be neglected since the primary WR star will 
dominate the EUV output. We feel that this ap- 
proach is a reasonable first order approximation 
to study the effect of WR formation through the 
binary channel on a stellar population. For more 
detailed studies the number of new parameters 
(mass ratio distribution, initial binary separation 
and more; see e.g. Dalton & Sarazin 1995, Van- 
beveren. Van Bever & De Bonder 1997) is large 
and several uncertainties associated with binary 
evolution (mass loss during Roche lobe overflow, 
mixing of secondary etc.; see e.g. Vanbeveren, 
Van Bever &: De Bonder 1997, Braun & Langer 
1995) should be taken into account. Population 
synthesis models which include binaries and are 
applicable to WR galaxies have also been pre- 
sented by Cerviho &: Mas-Hesse (1996), Cervifio, 
Mas-Hesse & Kunth (1996) and Vanbeveren, Van 
Bever & De Bonder (1997). First results from 



our work have been reported in Schaerer & Vacca 

(1996) . 

2.1.3. Emergent fluxes 

To describe the spectral evolution we rely on 
three different sets of theoretical models: 

1 ) For the main sequence evolution of massive 
stars we use the spectra produced by the com- 
bined stellar structure and atmosphere (CoStar) 
models of Schaerer et al. (1996a, b) and Schaerer 
& de Koter (1996). The models, which in- 
clude non-LTE effects, line blanketing, and stel- 
lar winds, cover the entire parameter space of 
O stars. As shown by Schaerer (1996b) and 
Schaerer & de Koter (1997) they represent sig- 
nificant improvements over previous calculations, 
particularly for the ionizing fluxes. In practice we 
use the model set given by Schaerer & de Koter 

(1997) for stars with initial masses Mini ^ 20 Mq. 
Additionally we restrict the use of these mod- 
els to the following T^s-logg domain: logg > 
2.2 and logg < 5.71 x logTeff - 21.95. The lat- 
ter restriction corresponds roughly to an 09/BO 
spectral type. This parameter range spans the 
domain of validity of the models and yields ion- 
izing fluxes which merge smoothly with those 
calculated from the Kurucz (1992) models de- 
scribed below. For solar metallicity we consis- 
tently use the corresponding CoStar models. At 
subsolar metallicities we adopt the low metal- 
licity {Z =0.004) calculations of Schaerer & de 
Koter. For Z =0.04 we use a new set of CoStar 
models (Schaerer 1996, unpublished). 

2) For stars less massive than 20 Mqwc use the 
line-blanketed plane-parallel LTE atmosphere mod- 
els of Kurucz (1992) calculated with a microtur- 
bulent velocity of Vturb = 2 km s~^. The Kurucz 
models cover a large range of metallicities, char- 
acterized by the iron abundance [Fe/H]. For each 
set of evolutionary tracks (characterized by a par- 
ticular value of the total metal abundance Z) we 
adopted the Kurucz model atmospheres with the 
corresponding value of [Fe/H]. The assignments 
for the specific metallicities considered in this pa- 
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per (dictated by the metallicities of the evolution- 
ary tracks) were as follows: we used the Kurucz 
models with [Fe/H] = -1.5 for Z = 0.001; [Fe/H] 
= -1 for Z = 0.004, [Fe/H] = -0.5 for Z = 0.008, 
[Fe/H] = for Z = 0.020, and [Fe/H] = -^0.3 for 
Z = 0.040. More accurate assignments are not 
necessary since the variations of the stellar spec- 
tra with metallicity are much less important than 
those due to changes of the evolutionary tracks. 

3) Stars in the WR phases are described by 
the spherically expanding non-LTE models of 
Schmutz, Leitherer & Gruenwald (1992). We 
identify the core radii and core temperatures of 
their models with the radius and Tefj from the 
wind-free hydrostatic stellar model (see Schaerer 
1996a for distinction and a discussion). Test 
calculations using the subphotospheric radii of 
Schaerer (1996a) yield essentially identical re- 
sults. The additional parameters which con- 
trol the emergent spectrum of WR stars are the 
mass loss rate M (taken from the stellar evolu- 
tion models) and the terminal velocity v^o of the 
wind. For WN and WC stars we use the fitting 
formulae from Schaerer (1996a, Eqs. 8 and 9), 
which account for the strong dependence of Voo 
on the WR subtype. For WO stars we adopt 
00=4000 km s~^, the average value of W03 and 
W04 stars (Kingsburgh & Barlow 1995, Kings- 
burgh, Barlow & Storey 1995). It should be 
noted that the atmosphere models of Schmutz, 
Leitherer & Gruenwald (1992) used here include 
only He and do not account for variations in 
the emergent spectra due to possible composi- 
tion or metallicity effects. These authors argue 
that uncertainties in the hydrodynamic struc- 
ture should have larger effects than the neglect 
of metals. For WN stars, pilot studies carried 
out by Schmutz (1991) indicate that the effects 
of line blanketing in WR atmospheres produce 
moderate changes in the predicted spectra, al- 
though Crowther et al. (1997) find a significant 
reduction of the ionizing flux for a line-blanketed 
WNL model. For WC/WO stars models which 
incorporate the abundant species of C and O (cf. 



Grafener et al. 1997) have just been developed. 
Once line-blanketed atmosphere models for WR 
stars have matured and become widely available, 
they will be incorporated in our synthesis rou- 
tines. 

In addition to the stellar continuum, we have 
included the nebular continuum spectrum in our 
models. As first shown by Huchra (1977) the 
nebular contribution becomes important when 
hot stars providing a large number of ionizing 
photons are present. To treat the nebular emis- 
sion in a simplified way we assumed the emitting 
gas has an electron temperature Tg = 10000 K, 
an electron density = 100 cm~^, and a helium 
abundance of 10 % by number relative to hydro- 
gen. The monochromatic luminosity of the gas, 
which is proportional to the number of Lyman 
continuum photons Qq, is given by 

Lx = ^^f,Qo, (2) 

where is the case B recombination coefficient 
for hydrogen, and is the fraction of ionizing 
photons absorbed by the gas. For the continous 
emission coefficient 7totai we used the atomic 
data from the tables of Aller (1984) and Fer- 
land (1980) for the wavelength regimes shortward 
and longward of 1 //m, respectively, to account 
for free-free and free-bound emission by hydrogen 
and neutral helium, as well as for the two-photon 
continuum of hydrogen^. 

2.2. Input parameters 

Evolutionary synthesis consists of calculating 
the properties of an entire stellar population as 
a function of time. In this work we concentrate 
on the time evolution of young populations expe- 
riencing star formation on timescales short com- 
pared to the evolutionary timescales of massive 
stars. As an idealized case we assume an instan- 
taneous burst occurring at time t = 0. The time 

■^We ignore two-photon continuum emission for A > 1 fj,m 
and assume 7hcI = 7h- 
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evolution of the burst is followed by isochrone 
synthesis based on a subroutine kindly provided 
by Georges Meynet. 

The parameters defining the burst population 
are the following: the initial metallicity Z, the 
initial mass function (IMF), and the total mass 
of stars formed in the burst. The last parameter 
serves only as a normalization constant and is 
not relevant for the quantities discussed in this 
work. 

For the IMF we adopt a power law (j)(m) = 
dN/dm oc between the upper and low cut- 
off masses, Mup and Miow respectively. The 
IMFs in massive-star forming regions in the Local 
Group generally have slopes (a) between 2 and 
2.8 (see e.g., Maeder & Conti 1994; Massey et al. 
1995b), in agreement with the Salpeter (1955) 
value of a = 2.35. We adopted the Salpeter IMF 
for our "standard" model. For the quantities and 
burst ages analyzed in this work the influence of 
the lower mass cut-off is negligible, and we use 
-^low = 0.8 Mq. a large number of observations 
indicate that the most massive stars have masses 
of the order of 100 M0or larger (e.g. Kudritzki 
1996). We adopt M^p = 120 M© for ah calcu- 
lations. For an accurate description of the evo- 
lutionary status of stars close to the lower mass 
limit for WR formation, Mwr, we use the tabu- 
lated values from Maeder & Meynet (1994). 

Results from calculations using variations of 
the "standard" input parameters are not pre- 
sented here, but are available in electronic format 
(see Appendix). 

2.2.1. WR emission lines 

To predict the broad stellar emission lines we 
have compiled average line fluxes for Of, WNE, 
WNL, WC, and WO stars. Our compilation, 
given in Tables |l] and ^, includes the strongest 
emission lines of H, He, C, and N in the UV and 
optical spectral ranges: He II A1640, the N iii/v 
A4640 blend, C iii/iv A4650, He II A4686, the 
He II A4861+H/3 blend (denoted as 4861), the 



He n+C IV A5411 blend, C III A5696, the C IV 
AA5808-5812 multiplet (denoted as 5808), and 
the He II A6562 -^Ha blend (denoted as 6560). 
All other optical emission lines from WR stars 
are expected to be much weaker, and therefore 
have not been included in our model calcula- 
tions. For comparison with previous models we 
also synthesize the two broad "WR bumps" (the 
N m/w+C iii/iv-hHe II blend at A ~ 4650, and 
C IV A5808). The cahbration from Smith (1991) 
of these features is reported in Table ^ for com- 
pleteness. 

For Of and WN stars we use the absolute flux 
in He II A4686 (col. 9 of Table ||) reference 
and compute line strength ratios for all other 
lines relative to this reference. For WC and WO 
stars we use C IV A5808 (col. 10 of Table P) as 
the line strength reference. All line fluxes are 
expressed in terms of these "reference lines". 

Of stars: The strongest emission lines in Of 
stars are He II 4686 and the N iii/v 4640 blend. 
The He II/N III ratio shows considerable varia- 
tions between individual Olf stars. We adopted 
a ratio of L(4640)/L(4686) = 0.1, which should 
strictly be considered a lower limit. For the He II 
4686 line we adopted an average luminosity de- 
termined from one 03If/WN6 and one 06Iaf star 
in the LMC (Crowther 1996, private communica- 
tion). This yields a value of 2.5 x 10'^^ erg s~^. 
This is a factor 2.6 lower than the line flux from 

Table 3: Average "WR bump" line luminosities 
for WN and WC stars. Logarithm of line lu- 
minosities in ergs~^. From Smith (1991, LMC 
stars) 



C iii/iv4650-hHe II4686 C IV5808 



WN7 36.5 35.1 

WC4 36.7 36.5 
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Table 1: Adopted line luminosities for broad stellar emission lines from WN stars. Mean absolute line 
luminosities are given for the reference line He II 4686 in erg s~^ for each subtype. Mean fluxes ratios 
with respect to this reference line are given for the other lines. Below each mean value are the standard 
deviation and the number of stars included in the mean. 



Subtype 




1640 
4686 


4640 
4686 


4861 
4686 


5411 
4686 


5808 
4686 


6560 
4686 


4686 


WNE 


Mean 


7.95 




0.129 


0.137 


0.074 


0.202 


5.2 X 10''^^ 




a 


2.47 




0.031 


0.022 


0.027 


0.130 


2.7 X 10^5 




N 


28 




7 


38 


5 


7 


26 


WNL 


Mean 


7.55 


0.244/0.616'' 


0.309 


0.107 


0.062 


1.596 


1.6 X 10''^^ 




a 


3.52 


0.151/0.519 


0.107 


0.022 


0.036 


1.511 


1.5 X 10^^ 




N 


17 


5/7 


11 


30 


10 


4 


19 



" The first set of values for N III 4640 were derived from LMC WNL stars and were used for models with 
Z < Zq; the second set of values were derived from Galactic WNL stars and were used for models with 
Z>Zq. 



the four 03f/WN stars in R136a recently ana- 
lyzed by de Koter, Heap & Hubeny (1997); for 
those stars L(4686) = 6.5 x 10^^ erg s"^ (de Koter 
1996, private communication). As can be seen 
the line luminosity of the R136a stars is com- 
patible with the average value derived for WNE 
stars. 

WN stars: Unfortunately, there is no single 
reference that provides intrinsic fluxes for all the 
lines in WN spectra that we wish to include in 
our models. We compiled flux measurements 
from the studies of Morris (1995) and Smith, 
Shara, & Moffat (1996, hereafter SSM96). Al- 
though Morris (1995) and SSM96 present results 
for many of the same stars, the two sets of mea- 
surements have only two lines in common: He II 
A4686 and He II A5411. We used the 5411/4686 
line ratio and the absolute value of the 4686 flux 
to check that the two sets of measurements were 
consistent. 



We began by dereddening the observed fluxes 
for the Galactic stars in both samples with the 
color excess values given by Morris et al. (1993), 
with the revisions given by Morris (1995). For 
three stars for which values of E(B — V) were not 
given we adopted the values given by Hamann et 
al. (1993). We used the reddening curve given by 
Seaton (1979) and Ry = 3.2. For the LMC stars 
we used the extinction curve given by Seaton 
(1979) for the Galactic foreground reddening and 
the average LMC extinction curve derived by 
Vacca (1992: sec also Morris et al. 1993) with 
Ry = 3.2 for the reddening within the LMC. We 
assumed a maximum foreground reddening value 
of E{B -V)fg = 0.03 mag. 

We adopted the subtype classifications given 
by SSM96. This rarely resulted in a change of 
more than one subtype with respect to older 
classifications (e.g., Conti & Massey 1989). We 
then divided the sample into WNE and WNL 
subclasses. The WNE subclass consists of stars 
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Tabic 2: Adopted line luminosities for broad stellar emission lines from WC/WO stars. Mean absolute 
line luminosities are given for the reference line C IV 5808 in erg s~^ for each subtype. Mean fluxes ratios 
with respect to this reference line are given for the other lines. Below each mean value are the standard 
deviation and the number of stars included in the mean. 



kJ LIU t V 




1640 
5808 


4650 
5808 


4686 
5808 


4861 
5808 


5411 
5808 


5696 
5808 


6560 
5808 


5808 




WO 

VV \J 


IVlcdli 




U.O-L 


u.oou 




n n9fi 

u.uzu 




U.U / o 


1.1 


X 


ir,36 




(T 


1 44 

J. 


n 18 


n 1 ^8 

U. -LOO 




n ni9 




Ofi5 
u.uuo 


n 9 
u.z 


s/ 

A 


1^36 

lU 




N 


2 


3 


3 




3 




2 


q 






WCA 

VV V^t: 




9 14 


1 71 




01 Q 


n 093 


n 091 


045 


3.0 


X 


iq36 




(7 


1 HQ 


\J.OO 




OOQ 

U.UUt/ 


OOQ 


U.U-LO 


U.U -L^ 


1 1 
1.1 


X 


1f>36 

lU 




N 


16 


23 




11 


14 


11 


19 


lc5 






VV 




^ 14 


9 




099 


n3fi 


U.U 1 1 


074 

U.U 1 ^ 


9.8 


X 


iq35 

lU 




(J 


4 4^ 


n 4Q 




01 ^ 

U.U -LO 


n ni8 


yj.yjoo 


01 1 

U.U-L ± 


2.3 


X 


1^35 

lU 




AT 

iv 


r 



y 




( 


Q 

o 


y 


Q 

o 


2 






WC6 


Mean 


7.60 


2.98 




0.047 


0.054 


0.184 


0.112 


8.9 


X 


1035 




a 


6.17 


0.52 




0.012 


0.006 


0.083 


0.037 


1.9 


X 


1Q35 




N 


6 


10 




9 


9 


10 


10 


3 






WC7 


Mean 


10.22 


3.21 




0.069 


0.076 


0.579 


0.177 


1.4 


X 


1q36 




a 


5.25 


0.49 




0.024 


0.027 


0.214 


0.028 


0.6 


X 


1q36 




N 


10 


15 




13 


14 


15 


14 


4 






WC8 


Mean 


13.69 


3.30 




0.096 


0.091 


1.968 


0.580 


3.5 


X 


10^5 




a 


14.58 


0.45 




0.037 


0.031 


0.649 


0.284 


3.1 


X 


1Q35 




N 


5 


7 




5 


6 


7 


8 


3 






WC9 


Mean 


4.47 


4.46 




0.294 


0.135 


3.705 


1.243 


2.0 


X 


10^5 




a 


1.25 


0.97 




0.095 


0.033 


0.514 


0.256 


2.0 


X 


1q35 




N 


3 


17 




17 


17 


18 


16 


2 
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with subtypes between WN2 and WN4, while the 
WNL subclass consists of stars with subtypes be- 
tween WN6 and WN9. This division agrees with 
the definitions and assignations given by Conti 
h Massey (1989). Consistent subclass identifi- 
cations could not be assigned to the WN5 sub- 
type and these stars were excluded from the aver- 
ages. We then computed the average dereddened 
5411/4686 ratios for LMC and Galactic stars for 
the two subclasses in each sample. It was found 
that the averages for the Morris sample and those 
for the SSM96 sample could be brought into 
agreement only if stars with low hydrogen abun- 
dances (WNb and WNo stars, as classified by 
SSM96) were excluded from the WNL averages. 
WR 43 was the only exception to this rule of ex- 
cluding objects classified by SSM96 as having low 
H abundance from the WNL averages; although 
SSM96 classifiy this object as a WN6o, Drissen et 
al. (1995; see also Dessart & Crowther 1997) find 
that it is composed of three WN6h stars. Cor- 
respondingly, we excluded those stars contain- 
ing a substantial amount of hydrogen (WNh or 
WN(h) stars, again as classified by SSM96) from 
the WNE averages. Once the mean line ratios 
were found to be consistent between the two sam- 
ples, we combined the samples. The flux ratios of 
He II A1640/He II A4686 and He II A5411/He II 
A4686 are taken from Morris (1995), while the 
other line strength ratios are taken from SSM96. 

As shown by SSM96 the strength of the N 
A4640 line depends on the initial abundance of 
the progenitor star. For solar metallicities and 
larger we adopt the Galactic average, while for 
Z < Zq we adopt the lower LMC value. For all 
of the other line ratios, we derived average val- 
ues by combining the LMC and Galactic stellar 
samples, after confirming that the means of the 
two samples for each line ratio were in agreement. 
We excluded from the averages any values which 
were found to be more than 3cj away from the 
mean; the averages were then re-computed with 
the smaller sample. This process resulted in the 
exclusion of only a few individual line flux values 



from the calculations of the various means. 

The absolute luminosity of He II for each 
subclass was derived from the Morris (1995) 
sample of LMC stars and those Galactic stars 
that are definite or probable members of clus- 
ter or associations with reliable distance esti- 
mates. We assumed a distance of 50.1 kpc 
(DM=18.50) for the LMC; distances for stars 
in Galactic clusters or associations were taken 
from Lundstrom & Stenholm (1984), with revi- 
sions given by Moffat (1983; for NGC 3603), Gar- 
many & Stencel (1992), Smith, Meynet, & Mer- 
milliod (1994), and Massey, Johnson, & DeGioia- 
Eastwood (1995a) and references therein. For 
WR 43/NGC 3603, we divided the total dered- 
dened flux measured by Morris (1993) into three 
components, with relative contributions in agree- 
ment with the line measurements of Drissen (1997, 
private communication). The resulting mean 
WN line luminosities and line ratios, as well as 
the standard deviations about the mean and the 
number of stars included in each mean, are pre- 
sented in Table (0). 

A comparison of the line fluxes determined 
for individual stars by Morris (1995) with those 
measured by Crowther (1997, private communi- 
cation) and Crowther & Dessart (1997) on sub- 
stantially higher resolution spectra reveals a sys- 
tematic difference of ~ 20 — 30 % (in the sense 
that the line fluxes from Morris 1995 are larger 
than those of Crowther & Dessart 1997). It is 
not completely clear what the cause of this dis- 
crepancy might be, but if we assume that the 
higher resolution spectra yield more accurate line 
fluxes, then we conclude that our adopted line 
luminosities may be ~ 30 % too large. In ad- 
dition, the absolute luminosity of He II A4686 
exhibits substantial variations within the WNL 
subclass, and a mean luminosity for this line is 
rather uncertain. The situation is illustrated in 
Fig. "Kl, in which we plot the line luminosities 



from Crowther (1997) and Crowther & Dessart 
(1997) against the bolometric luminosities de- 
rived by these authors from their spectral analy- 
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ses. The large variation in line luminosity within 
the WNL subclass is easily seen. The figure also 
reveals that there may be a difi'erence in the line 
luminosity between stars with L < IO^Lq and 
those with L > IO^Lq. The lower luminosity 
stars (in the combined Galactic plus LMC sam- 
ple) have a mean He II 4686 line luminosity of 
5.6 X 10^^ erg s~^, while the higher luminosity 
stars have a mean line luminosity of 3.1 x 10^^ 
erg s~^. However, the uncertainties in the mean 
values for these two luminosity classes are also 
very large, and we have chosen to use a constant 
line luminosity for all bolometric luminosities of 
WNL stars in our models. The entire WNL sam- 
ple from Crowther (1997) has a mean line lumi- 
nosity of 1.1 X 10^^ erg s~^, which agrees well with 
our adopted value after the systematic differences 
between the Crowther (1997) and Morris (1995) 
line flux measurements are taken into account. 

Despite the uncertainties in the mean lumi- 
nosity of the He II 4686 line, our estimate of 
1.6 X 10'^^ erg s~^ is in excellent agreement with 
the value of 1.7 x 10^^ erg s~^ determined by 
Vacca (1992; see also VC92). Vacca (1992) 
demonstrated that this value, when combined 
with spatially- integrated line flux measurements, 
yields an estimate for the number of WNL stars 
within the 30 Doradus complex that is in good 
agreement with the number derived by directly 
counting objects with spectroscopic identiflca- 
tions. Furthermore, Terlevich et al. (1996) 
used this mean luminosity value to determine 
the number of WN stars in NGC 604 from their 
spatially-integrated spectra. They found that 
their estimate was again in excellent agreement 
with the number of WR stars found by Drissen, 
MoflFat, & Shara (1993) using narrow-band HST 
images. In addition, Crowther &: Dessart (1997) 
flnd that this value of the mean line luminosity 
yields a number for the WNL stars within the 
NGC 3603 complex that is in good agreement 
with the observed number (Drissen et al. 1995). 
Unfortunately, additional tests of our value are 
currently not available. 



To reduce the uncertainties of the WR line 
strengths used in the present models it would be 
highly desirable to understand the observed large 
variations of the line luminosities (Fig. 1) and 
their dependence on stellar parameters. If corre- 
lations of the line strengths with physical param- 
eters can be established, these can be used in fu- 
ture models to provide a more accurate picture of 
the WR populations. It is also necessary to clar- 
ify the relationship between WR classiflcations 
and stellar models (cf. Smith & Maeder, 1997). 
For example, the WNE and WNL classiflcations 
we have adopted are deflned according to purely 
phenomenological criteria and do not necessarily 
correspond directly to the deflnitions employed in 
the stellar evolutionary models. Our additional 
separation of the stellar sample into objects with 
and without hydrogen represents a compromise 
between the observational definitions of WNE 
and WNL stars and those employed by the evo- 
lutionary models. With these criteria we have 
attempted to link the observations and models as 
closely as possible. The adopted approach should 
yield average line luminosities and ratios that are 
representative of the properties of the classes of 
hydrogen-free and hydrogen-rich WN stars. 

WC and WO stars: In view of the substan- 
tial variations in the frequency distributions of 
WC/WO subtypes with environment (cf. Smith 
& Maeder 1991), as well as the strong variations 
in spectral signatures with subtype, we chose to 
treat the various WC subtypes separately. The 
assignment of the WC/WO subtype from the 
evolutionary model is discussed below (§^^). 

The task of compiling dereddened line fiuxes 
was substantially easier in the case of WC and 
WO stars. Average line fiuxes of WC stars have 
been derived by Smith, Shara, & Moffat (1990a, 
b, hereafter SSM90a, SSM90b) and recently by 
Brownsberger (1995), who provided an extensive 
list of dereddened line fiuxes for nearly every 
known WC star in the Galaxy and the LMC. 
Their results agree within their quoted uncer- 
tainties. We adopted the dereddened line fiuxes 
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presented by Brownsberger (1995) and computed 
average values of the line flux ratios (relative to 
the strength of C IV A5808) function of 

WC subtype; the subtypes were assigned accord- 
ing to the classification scheme of SSM90a and 
SSM90b. (The method used by Brownsberger to 
compute the reddening for each star and dered- 
den the line fluxes is identical to the method we 
adopted for WN stars; see above.) The average 
luminosity in C IV A5808 for each subtype was 
determined from WC stars in the LMC and those 
Galactic stars located in clusters or OB associa- 
tions for which reliable distances could be deter- 
mined. Distances were taken from Lundstrom & 
Stenholm (1984), with revisions given by Gar- 
many &; Stencel (1992), Smith, Meynet, &: Mer- 
milliod (1994), and Massey, Johnson, & DeGioia- 
Eastwood (1995a) and references therein. Unfor- 
tunately it is not possible to check for consis- 
tency in the C IV line luminosity between LMC 
and Galactic WC stars; the LMC contains only 
one WC star with a subtype later than WC4, 
and there are no WC4 stars in the Galaxy with 
well determined distances. For the WC4 subtype 
we adopted the average C IV A5808 luminosity 
of the LMC stars. This value is 0.5 dex larger 
than that found for WC5-7 stars in the Galaxy, 
which agrees with the findings of SSM90b. The 
mean line ratios were derived by combining the 
LMC and Galactic WC4 samples, after confirm- 
ing that their individual mean values were com- 
pletely consistent with one another. For all other 
types, we adopted the mean line luminosity and 
line ratios derived from the Galactic stars. The 
average line luminosities and ratios, as well as 
the standard deviations about the means and the 
number of stars included in the estimation of the 
means, are given for each subtype in Table |2[ The 
averages were computed iteratively, as described 
above for the WN stars, by excluding those indi- 
vidual values that were found to be more than 3 
a from the means. This resulted in the exclusion 
of only a few line flux values from the calculations 
of the various means. 



The emission of WC stars in the 4650 bump 
requires an additional comment. The broad 
C III/C IV A4650 blend also includes the con- 
tribution from He II A4686, which usually can 
not be separated. From SSM90b (their Table 2) 
we find that for WC5-9 stars He II 4686 con- 
tributes on the average ~ 12 % to the 4650 blend; 
SSM90a estimate He II 4686 contributes 8 to 30 
% of the blend for WC4. We include the strength 
of the entire 4650 blend in our models and do not 
treat the He II 4686 emission from WC stars sep- 
arately. 

For WO stars we have derived average line 
luminosities from the dereddened line flux mea- 
surements of three W03 and W04 stars given 
by Kingsburgh & Barlow (1995) and Kingsburgh, 
Barlow & Storey (1995). While the line luminosi- 
ties of these objects (all located in low metallicity 
metallicity environments: LMC, SMC, IC1613) 
exhibit relatively little scatter, they are system- 
atically stronger than the line luminosities for 
the two known galactic WOl and W02 stars. 
This difference may be attributed to differences 
in their evolution (Kingsburgh, Barlow & Storey, 
1995). Since WO stars are expected to occur 
mainly at low metallicities (cf. § |3.2D it seems jus- 
tified to adopt the properties of only W03 /W04 
stars. We used the fractional line contributions 
provided by Kingsburgh, Barlow &: Storey (1995) 
and Kingsburgh & Barlow (1995) to decompose 
the total fiux in the emission blend at 4650 A 
into the flux from C IV A4658 and that from 
the component due to C IV A4685, and He II 
A 4686. (The latter two lines would not usu- 
ally be separable in the spectra of W-R galaxies.) 
From test calculations we flnd that other emis- 
sion lines which are prominent only in the rare 
WO stars are predicted to be very weak and are 
therefore not included in our models. For com- 
pleteness we list their average line luminosities, 
which are: L(0 IV A3400) = 1.9 x lO^^ erg s-^ 
L(0 VI A3811) = 1.5 X 10^^ erg s'^, and L(0 V 
A5590) = 1.5 X 10^^ erg s~^. Note that we have 
adopted a distance of 58.9 kpc for the SMC and 
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661 kpc for IC 1613. 

2.3. Synthesized quantities 

The most important synthesized quantities 
can be classified in three groups: population 
statistics, nebular quantities, and WR emission 
lines. 

2.3.1. Population statistics 

For comparisons with star counts it is use- 
ful to derive the numbers of stars with different 
spectral types. Of particular interest for young 
bursts is the number of O stars, the total num- 
ber of WR stars and their distribution among 
the WNL, WNE, and WC/WO subclasses. WR 
stars and the various WR subclasses are defined 
by their surface abundances and their effective 
temperature from the evolutionary models (see 
e.g. Maeder 1991b). Hydrogen-free WN mod- 
els are considered to be WNEs, while those con- 
taining hydrogen are treated as WNLs. We also 
estimate the number of stars belonging to the 
various WC subtypes (WC4 - WC9) and WO 
subclasses (cf. below); the frequency distribu- 
tion among these subtypes depends strongly on 
metallicity (e.g. Smith & Maeder, 1991). We as- 
sign WC and WO subtypes according to their 
surface abundance, i.e. by the (C-|-0)/He num- 
ber ratio. Note that in some cases the WR 
phase can be entered during the core H-burning 
phase (cf. Meynet et al. 1994). We define 0-type 

WC/WO subtypes are assigned as follows: WC9: 
(C+0)/He < 0.08, WC8: (C+0)/He G (0.08,0.15], WC7: 
(C+0)/He e (0.15,0.25], WC6: (C+0)/He G (0.25,0.35], 
WC5: (C+0)/He € (0.35,0.43], WC4: (C+0)/He G 
(0.43,1.], WO; (C+0)/He > 1. Here the values of Smith & 
Maeder (1991) have been adjusted slightly to account for 
the recent C, O, and He determinations in WC4 stars by 
Grafener et al. (1997). The correspondence between spec- 
tral types and C, O abundance is not completely agreed 
upon, however (cf. e.g. Eenens & Williams 1992 and 
Koesterke & Hamann 1995). More detailed analyses, cov- 
ering a large range of WC subtypes and including oxygen 
abundance determinations, will be necessary to assess the 
validity of the correlation of (C-|-0) /He with WC subtype. 



stars to be those H-burning objects with effective 
temperatures Tefj > 33000 K. 

We also estimate the number of Of stars and 
include the contribution of their stellar emission 
lines to the observed emission features. We as- 
sumed that all O stars with surface gravities be- 
low the average value found for Ola stars can be 
counted as Of stars. We used the recent cali- 
bration of Vacca et al. (1996) to derive log 5 as 
a function of T^s for the O la luminosity class. 
For consistency we adopt the calibration given 
for the gravity 5evol determined from evolution- 
ary models (following Vacca et al. 1996). 

2.3.2. Nebular quantities 

We synthesize the following H and He nebu- 
lar recombination lines under the assumption of 
Case B recombination: Hq, H/5, He I A 4471, 
and He II A 4686. The adopted electron temper- 
ature and density are the same as given above. 
The emissivities and recombination coefficients 
are taken from Osterbrock (1989). The lumi- 
nosity emitted in the emission lines is calculated 
from the number of Lyman, He I, and He II pho- 
tons (Qo) Qi) and Q2 respectively) using the re- 
lations 

L(Ha) = 1.36 X lO^^ViQo (3) 
L(H/3) = 4.76 X 10-^^7 Qo (4) 
L(4471) = 2.22 X 10-^^7 <9i (5) 
L(4686) = 1.02 X 10"^V7Q2 (6) 

where L is given in [ergs""*^] and Qi in [s~"^]. /-y is 
the fraction of ionizing photons absorbed by the 
gas. A fraction (1 — f^) is assumed to escape the 
system. Our "standard" model assumes an ion- 
ization bounded nebula with = 1. It should be 
noted that the strength of He I A 4471 predicted 
from Eq. ^ above will be incorrect when there are 
enough photons above 24.6 eV to fully ionize he- 
lium in the nebula; this situation corresponds ap- 
proximately to equivalent effective temperatures 
of Teff > 40000 K. 
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2.3.3. WR emission lines 

Several observational signatures of WR stars 
in the integrated spectrum of a starburst re- 
gion are computed according to the prescriptions 



given in Sect. p. 2.1 



3. Wolf-Rayet and O star populations in 
young starburst 

We now discuss the model predictions regard- 
ing the massive star populations. We present the 
time evolution of the WR and O star populations 
for an instantaneous burst at different metallici- 
ties. We explicitly consider the various WR sub- 
types (primarily WNL and WC/WO) and esti- 
mate the fraction of O stars which exhibit emis- 
sion lines (Of stars). The results presented here 
are essential for understanding the behaviour of 
the synthesized spectral lines and other observ- 
able features which will be discussed in § ^ and 



In Figure we show the time evolution of the 
stellar number ratios WR/(WR+0), WNL/(WR-F0), 
and WC/(WR+0), as well as the Of/0 ratio 
for the five metallicities between Z=0.001 and 
0.040. (Note that here WC stands for WC and 
WO stars, irrespective of the distinction made 
between these types in the calculation of the line 
strengths.) The values were calculated for an 
instantaneous burst (occurring at t = 0) and 
a Salpeter IMF. Variations in the relative WR 
and O star populations with other parameters 
(IMF slope, upper mass limit, star formation rate 
and different evolutionary tracks) have been dis- 
cussed in detail by Meynet (1995), and therefore 
will not be repeated here. Nevertheless, we must 
comment on a small discrepancy between our re- 
sults for WR populations and those presented by 
Meynet (1995). The differences in the two sets 
of results can be seen by comparing, for exam- 
ple, our predictions for the WC population at 
Z=0.004 (Fig. ^) with those presented in Fig. 3 
from Meynet (1995). As Fig. shows, we pre- 
dict a strong decrease of the WC population after 



~ 4 Myr, while in models of Meynet (1995) the 
WC population remains dominant until the end 
of the WR-rich phase (t ~ 5.2 Myr). The ori- 
gin of this discrepancy resides in the use of dif- 
ferent techniques to interpolate between the two 
(discrete) evolutionary tracks above and below 
the mass limit of WR formation (in the case of 
Z=0.004, between the 60 and 40 Mq track). This 
discrepancy occurs therefore only for ages larger 
than the lifetime of the least massive star whose 
initial mass is larger the WR mass limit (60 Mq 
in the present case). While Meynet (1995) per- 
forms a linear interpolation on the lifetimes in the 
various WR phases (WNL, WNE, and WC) as a 
function of initial mass between those two tracks, 
our classification of WR subtype is based on an 
interpolation of the surface abundances on the 
isochrone. Meynet's method predicts the exis- 
tence of WC stars until the end of the WR phase, 
while in our case models close to the WR mass 
limit have surface abundances typical of the less 
evolved WNL stars, since mass loss is not suf- 
ficient to expose more processed matter. This 
numerical difference explains the shorter dura- 
tion of the WC phase (and hence correspondingly 
longer WNL phase) in all our models compared 
to the results of Meynet (1995). For the same 
reason the duration of the WNE phase in our 
models is significantly shorter. Although both 
the exact behaviour of stars in this mass range 
and the value of Mwr are somewhat uncertain, 
we feel that the present approach yields a more 
physical description of the transition between the 
least massive WR stars and the lower mass stars. 
The following important points are illustrated 



in Figure 



The fraction of Of/0 stars (classified by 
the criterion given in §^]^) is typically of 
the order of 0.1 to 0.3 for burst ages up 
to 2-3 Myr. The predicted changes with Z 
are a result of the differences in the main 
sequence as a function of Z. Due to the 
blueward shift of both the zero-age main se- 
quence and the terminal-age main sequence 
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Table 4: Maximum WR/0 ratios 



Z WNL/0 WR/0 

a = 1. a = 2. a = 2.35 a = 1. a = 2. a = 2.35 



0.040 


1.50 


0.54 


0.36 


1.50 


0.54 


0.36 


0.020 


0.51 


0.17 


0.11 


0.97 


0.32 


0.22 


0.008 


0.23 


0.12 


0.10 


0.40 


0.14 


0.11 


0.004 


0.13 


0.058 


0.044 


0.29 


0.091 


0.059 


0.001 


0.12 


0.033 


0.020 


0.15 


0.041 


0.025 



(t = 3 Myr) 
(t < 4 Myr) 



(cf. Meynet et al. 1994), the fraction of its 
main sequence-lifetime a star with an given 
initial mass spends as an O star increases 
with lower Z . This leads to the general de- 
crease of Of/0 to low Z . 

• The maximum of the WR/0 ratio (and 
therefore also the WNL/0 and WC/0 ra- 
tios) and the duration of the WR-rich phase 
decrease rapidly with metallicity (cf. Maeder 
& Meynet, 1994; Meynet 1995). The max- 
imum predicted values are given in Table 
|. Note that for Z < Zq our WNL/0 val- 
ues are slightly larger than those given by 
Meynet (1995, his Fig. 6) although both 
results are based on the same evolution- 
ary tracks. This difference is essentially 
due to the overestimate by Meynet of the 
WC population, as discussed above. For 
Z > 0.020, the WR/0 star number ratio 
formally becomes infinite at ages t > 4.5 
Myr because the O star population (defined 
to be those objects with T^s > 33000 K) 
disappears completely while WR stars are 
still present. Therefore, for these metal- 
licities in Table ^ we give values for the 
WR/0 ratio at specific times t < 4.5 Myr 
(cf. Meynet 1995). 

• In an instantaneous burst the WR-rich pe- 



riod (WR/0 > 0) is characterized by three 
distinct phases: 1 ) The early WNL phase 
at the beginning of the WR-rich phase oc- 
curs before the most massive stars evolve 
to WC stars. 2) The WC-rich or WC- 
dominated phase with coexisting WNL 
stars, where the precise mixture of differ- 
ent WR subtypes depends critically on age 
and metallicity. 3) The final WNL phase 
with WR stars which have evolved from ob- 
jects with initial masses close to the lower 
mass limit of WR formation Mwr. Phases 
occurring before and after the WR-rich pe- 
riod have been discussed by Meynet (1995), 
to which we refer the reader for more de- 
tails. 

3.1. Frequencies of WR-rich starbursts 

An interesting result from the work of Meynet 
(1995) is the predicted frequency of WR-rich ob- 
jects in young starbursts. In Table ^ we present 
the revised frequencies predicted from our mod- 
els. For an instantaneous burst scenario one can 
easily determine the time during which the popu- 
lation contains only O stars ("Ooniy")) both WR 
and O stars ("WR-hO"), or WR stars and no O 
stars (i.e. B instead of O stars: "WRoniy")- The 
values given in Table || (cols. 2-4) are the fraction 
of time spent in these phases with respect to the 
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Table 5: Statistics of H II and WR galaxies. The 
metallicity is given in col. 1. Cols. 2-6 indicate 
the predicted fraction of H II regions contain- 
ing: only O stars (col. 2), WR and O stars (3), 
WR but no O stars (4) . The fraction of H II re- 
gions dominated (in number) by WNL stars, and 
WC stars is given in Cols. 5 and 6. All numbers 
are relative fractions of the total number of pre- 
dicted H II regions defined as clusters harbouring 
O and/or WR stars. Col. 7 indicates the fraction 
of WR-rich H II regions dominated by WC stars 
(in number). 



z 


Oonly 


WR+0 


WRonly 


WNL 


WC 


WC 
totalWR 


0.040 


0.20 


0.26 


0.54 


0.48 


0.32 


0.40 


0.020 


0.28 


0.36 


0.36 


0.50 


0.22 


0.31 


0.008 


0.40 


0.60 




0.43 


0.17 


0.28 


0.004 


0.60 


0.40 




0.27 


0.13 


0.33 


0.001 


0.82 


0.18 




0.18 







total time when O and/or WR stars are present. 
Columns 5 and 6 give the fraction of time during 
which the WR population is expected to be domi- 
nated (in number) by stars of WNL and WC sub- 
types respectively. The last column (col. 7) gives 
the fraction of WR-rich starbursts that are ex- 
pected to be dominated by WC stars (i.e. the val- 
ues in column 7 are derived by dividing the values 
in column 6 by the sum of the values in columns 
3 and 4). Note that the values in cols. 5-7 differ 
from those of Meynet (1995, Table 1, high mass 
loss models) whose results overestimate the rela- 
tive duration of the WC and WNE phases with 
respect to the WNL phase as discussed above^ 

For an instantaneous burst of star formation, 
the phases during which O and/or WR stars are 
present in the population should lead to a strong 
emission line spectrum. Therefore, the fractions 
given in Table |5| should represent the fraction of 
H II regions (extragalactic H II regions or alike 
objects) which should be dominated by the var- 
ious stellar types. Although comparisons with 
observations are difficult due to problems regard- 
ing the statistical completeness of a given sample, 
estimates derived from low metallicity samples 
yield encouraging results regarding both the fre- 
quency of the WR phase and the occurrence of 
WC stars (Schaerer & Vacca 1996). Systematic 
studies of the WR and O content in statistically 
complete samples should be possible in the near 
future. 

3.2. WC/WO populations in starbursts 

From both observational and theoretical stud- 
ies it is well known that the subtype distribution 
of WC stars changes dramatically with metallic- 
ity (cf. Smith & Maeder, 1991). In particular, 
late type WC stars are found only in high metal- 
licity environments (e.g., Conti & Vacca 1991; 
Philips & Conti 1992), while all WC/WO stars 



^In Table 1 of Meynet (1995), columns 1 and 2 for the high 
mass loss models at Z=0.004 are misprints. All other 
values in common with Meynet are in agreement. 
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are found to be WC4 or WO types in low metal- 
licity regions. 

WC stars are expected to evolve towards ear- 
lier subtypes during their lifetime, an evolution 
that reflects the change in surface abundance as 
larger amounts of processed matter are brought 
to their surface (Maeder 1991b). Therefore in 
a burst at a given metallicity, the distribution 
of WC /WO subtypes is expected to evolve with 
time. To quantify this behaviour we have cal- 
culated the IMF-weighted mean WC subtype 
(henceforth called IMF-averaged subtype), which 
reflects the average WC subtype of a mixed popu- 
lation of WC stars. We assigned WO stars to the 
WC3 subtype. The time evolution of the IMF- 
averaged WC subtype at different metallicities is 
plotted in Figure |A3| . The horizontal lines on 
the left indicate the average value for regions of 
constant star formation. 



Figure A3 nicely illustrates the qualitative 
changes of the WC population not only with 
metallicity but also with time. As expected this 
figure shows the predominance of WO and early 
type WC stars at low metallicity. Similarly, at a 
given Z the average WC subtype decreases with 
time due to the increase of the surface (C+0) /He 
abundance ratio along each evolutionary track. 
Late WC stars are therefore expected only at 
early times during the WC-rich phase (typically 
at t ~ 3-4 Myr), or during later times in a high 
metallicity region. 

4. Evolution of the ionizing spectrum 

In this Section we discuss the behaviour of 
quantities directly related to the ionizing flux of 
the stellar population. 

4.1. Ionizing fluxes 



Figure A4 presents the predicted time evolu- 
tion of the photon luminosity in the ionizing con- 
tinua of hydrogen (Qo); He I (Qi), and He II 
(Q2)- Note that all values are normalised to a 
population of a total mass of 1 Mq assuming a 



Salpeter IMF from 0.8 to 120 Mq. However, the 
qualitative properties discussed below are appli- 
cable to most IMFs. 

The metallicity dependence of the ionizing 
continuum of an integrated population is deter- 
mined primarily by variations in the evolution- 
ary tracks; changes due to blanketing effects in 
the atmosphere models are relatively minor (e.g. 
Cervifio & Mas-Hesse 1994). The behaviour of 
the Lyman continuum flux is controlled essen- 
tially by the main sequence evolution. In WR- 
rich phases the maximum contribution to the Ly- 
man continuum provided by these stars is typi- 
cally ~ 30 % (cf. also below). Main sequence 
stars also govern the behaviour of Qi, although 
the picture is somewhat more complicated by the 
fact that WR stars contribute a non-negligible 
fraction of the He I ionizing flux. In very young 
bursts (< 3 Myr) the He II ionizing continuum 
is dominated by O stars. The small differences 
with metallicity in this case are due to changes 
in the stellar wind properties, which dominate 
the opacity in the He II continuum (see Schaerer 
& de Koter, 1997). During the WR-rich phases 
Q2 is clearly dominated by the contribution of 
WR stars. The time and metallicity dependence 
of Q2 is therefore directly related to the dura- 
tion of this phase and the total number of WR 
stars present (the "strength" of the WR phase 
illustrated in Fig. [A2|) . 

What implications do the recent stellar atmo- 
sphere models have for the predictions of the ion- 
izing flux of a stellar population? As discussed 
in detail by Schaerer & de Koter (1997) the in- 
clusion of non-LTE effects, line blanketing, and 
stellar winds in the CoStar models leads to sig- 
nificant changes in the ionizing fluxes of O stars. 
The predictions of these models appear to be sup- 
ported by observations of H II regions (Stasihska 
& Schaerer, 1997). Figure illustrates the dif- 
ferent predictions for the ionizing fluxes for an 
instantaneous burst at solar metallicity from the 
various atmosphere models for O and WR stars. 
The solid line shows our standard model, using 
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the CoStar models for O stars and the fluxes 
from Schmutz, Leitherer & Gruenwald (1992) for 
WR stars. The dashed line was derived using the 
Kurucz (1992) plane parallel, LTE model atmo- 
spheres models instead of the CoStar models for 
O stars. Adopting Kurucz model atmospheres 
for both O and WR stars (using the wind cor- 
rected effective temperature from the evolution- 
ary models) yields the dotted line. 

Figure ^ (cf. dotted and dashed) illustrates 
the importance of using appropriate atmosphere 
models for WR stars to describe the ionizing 
fluxes in WR-rich phases (between ~ 3 and 6 
Myr). The role of WR stars in this context is 
discussed in detail by Schmutz, Leitherer &; Gru- 
enwald and Garcia- Vargas (1996). As can be seen 
from Fig. |A5| (comparing dashed and solid lines) 
the inclusion of CoStar models for O stars leads 
to the following changes with respect to synthe- 
sis models relying on older (and less complete) 
stellar atmospheres: In young bursts Q2 in- 
creases by at least 4 orders of magnitudes. The 
absolute value of the ionizing flux in the He II 
continuum can therefore reach values similar to 
those attained during the WR-rich phase. 2) In 
the first few million years of a burst the He I ion- 
izing flux increases by typically 60 %. 3) Small 
differences in the Lyman continuum flux are due 
mostly to larger values of Qq for supergiants with 
M = 20 to 40 Mq. 

4.1.1. The ionizing flux of WR stars in star- 
bursts 

As mentioned above the contribution of WR 
stars to the ionizing spectrum of a stellar pop- 
ulation is not negligible, especially at the short- 
est wavelengths. In this context, it is interesting 
to examine the average Lyman continuum lumi- 
nosity emitted per WR star, Qo(WR)/iV(WR), 
where (5o(WR) is the total Lyman luminosity 
emitted from all WR stars and A^(WR) is the 
number of WR stars present. This quantity is 
plotted in Fig. where we show the temporal 
evolution of Qo(WR)/iV(WR) for all metalUci- 



ties. 

For each metallicity the average Qq emitted 
per WR star generally decreases with time as ex- 
pected from the decrease of their average lumi- 
nosity with progressing age. The increase of the 
maximum average Lyman luminosity (and also 
of the time-averaged Qo(WR)/iV(WR)) with de- 
creasing Z simply reflects the higher average lu- 
minosity of WR stars at low Z, or equivalently, 
the larger value of Mwr - the minimum mass 
of WR formation - at lower Z (cf. § ^^). It is 
interesting to note that the average value of Qo 
per WR star shows variations over typically more 
than one order of magnitude during the WR rich 
phase. As will be discussed below (§ ^ ) this 
variation affects estimates of WR/0 ratios de- 
rived from He II A4686 and H/3 line strengths; in 
the past such estimates were calculated assuming 
a single dominant WR subtype and a constant 
value of Qo per WR star (e.g. Qq = 1.7 x 10"^^ 
s-i for WNL stars, VC92). 

5. Evolution of nebular lines of hydrogen 
and helium 

5.1. H/3 

Hydrogen recombination lines, particularly H/3 
serve as a fundamental reference for measure- 
ments of other emission lines. In addition, under 
the assumption of an ionization bounded neb- 
ula, their luminosities provide an estimate of the 
ionizing flux present and their equivalent widths 
can be used as an accurate indicator of the age 
of the population in the case of an instantaneous 
burst of star formation (a coeval population), as 
first pointed out by Copetti, Pastoriza & Dottori 
(1986). A critical discussion of this and other age 
indicators for young starbursts can be found in 
Stasinska & Leitherer (1996). 

The temporal evolution of the H/3 equivalent 
width 14/^ (H/3) for all metallicities is shown in Fig. 
A7 (left panel). The right panel shows VF(H/3) 
only during WR rich phases. The results pre- 
sented in these figures can be combined with re- 
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suits presented below (cf. Figs. [AI4 |A15| ) to de- 
termine the WR and O star populations in a re- 
gion from observational quantities. 

5.2. He II A4686 



Figure AS shows the predicted nebular He II 
A4686 emission for all metallicities; this figure 
was presented in Schaerer (1996c). Here we 
briefly reiterate the most important points demon- 
strated in this figure. 

For young bursts dominated by O stars (t ^ 3 
Myr) typical values of /(He Il)//(H/3) are be- 
tween 5 X 10""^ and 2 x 10~^. As mentioned 
above (cf. § due to the strong He II continuum 
flux obtained from the most recent and physi- 
cally complete O star models, these values are 
~ 4 order of magnitudes larger than predictions 
from other current synthesis models. However, 
the /(He II)//(H/?) values are still too low to be 
detectable in nebulae around O stars or in young 
H II regions. 

During the subsequent WR rich phase /(He II) 
/ I(H/?) values of IQ-^ to lO^^ are obtained. 
In order to determine which sources contribute 
to the high value of /(He Il)//(H/3) we have 
plotted the time evolution of 4686/H/3 during 
the WR phase for two selected metallicities in 
Fig. (upper left panel: Z=0.004, upper right: 
Z=Omi) and the corresponding WR/(WR-hO), 
WNL/(WR+0), and WC/(WR+0) number ra- 
tios (lower left panel: Z=0.004, lower right: 
Z=0.001). This Figure indicates that large neb- 
ular He Il/H/3 values are obtained during the 
WC/WO phase as well as during a short time 
before and after the WC phase (cf. Fig. and 
Schaerer 1996c). The high excitation is provided 
by WC/WO stars and WN stars that are about 
to evolve to the former phase. In addition, the 
relative contributions of WN and WC/WO stars 
are a strong function of both age and metallic- 
ity. This is illustrated by presenting the value of 
/(He Il)//(H/3) obtained when the contribution 
of WC/WO stars to the ionizing flux is excluded. 
For Z=0.004 almost no change in the line inten- 



sity ratio is seen because the nebular He II emis- 
sion is due primarily to hot WN stars with low 
hydrogen surface abundances (H mass fraction 
< 0.1). For Z=0.001 the high /(He Il)//(H/3) 
ratio is due mainly to highly evolved WC/WO 
stars (in this case primarily the WO subtypes, 
cf. Fig. which have a large Q2/Q0) despite 
their small numbers relative to WN stars. In con- 
clusion, both WC/WO and hot WN stars pro- 
duce the nebular He II emission. Their relative 
contributions, however, depend very sensitively 
on metallicity and (especially) age. 

The measurement of the observed flux in a 
nebular He II A 4686 line from a starburst region 
is complicated by the possible presence of a broad 
stellar component. In some cases, the narrow 
nebular component will be superposed on a broad 
feature (cf. e.g. VC92). However, in many in- 
stances, the broad feature may completely domi- 
nate the emission, which makes the identification 
and measurement of a nebular component nearly 
impossible with low to medium resolution spec- 
troscopy. The interplay between the two compo- 
nents will be discussed in S 6.2.3. 



6. Evolution of broad stellar emission lines 

In this Section we discuss the behaviour of 
the broad emission lines included in our calcu- 
lations (see § 2.2.1| ). Two different types of illus- 
trations are useful for comparisons with observa- 
tional data: 1 ) Plots of the equivalent widths of a 
feature as a function of age, and 2) Line intensity 
ratios with respect to H/3 function of the H/3 
equivalent width, which serves as an age indica- 
tor. Depending on the spectroscopic observations 
(type of object, spatial resolution, quality of ob- 
servations etc.), and several other factors one or 
the other approach may be favoured. Figures of 
the first type depend only on stellar quantities, 
since in most cases nebular continuum emission 
is fairly small. Equivalent widths may, however, 
be reduced by an underlying population. Inten- 
sity ratios of the WR lines with respect to H/3, 
for example, have the advantage of being rela- 
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tively independent of such dilution effects. On 
the other hand comparisons of this type rely on 
the assumption of an ionization bounded nebula, 
and may be affected by the relative spatial dis- 
tribution of the stellar and nebular emission, and 
the extent of the emission compared to the size 
of the aperture used for the spectroscopic obser- 
vations. 

6.1. UV lines: He II A1640 



total "bump" necessary in many cases, although 
working with the individual WR lines would be 
preferred. 



Figure AlC shows the temporal evolution of 
the He II A1640 equivalent width for all metal- 
licities. As with all quantities related to WR 
stars, the substantial variation with metallicity 
is due to the strong dependence of WR evolution 
on Z (cf. Maeder & Meynet 1994). For subso- 
lar metallicities 1^(1640) is always predicted to 
be smaller than ~ 4 A. For these metallicities, 
the time dependence of TV (1640) follows the total 
WR/ (WR-l-0) ratio quite closely, as discussed in 
§ ^.3.3| . Since for Z > 0.02 WN stars dominate 
the 1640 emission (late WC stars have a lower 
1640 emission), stronger variations are found be- 
tween the WC and WN dominated phases. 

6.2. Optical lines 

6.2.1. The traditional WR bumps at low spec- 
tral resolution 



In Figure All we plot the evolution of the 
equivalent width of the two WR emission blends 
prominent in the optical wavelength regime; these 
are the strongest WR features found in the spec- 
tra of WR galaxies. Adopting the line fluxes 
from Tables |l] and ^ we have calculated the 4660 
A bump (solid line) from the sum of N lii/v A 
4640, C iii/iv A 4650, and He II A 4686; the 
5808 A emission blend (long-dashed) is due to the 
C IV multiplet. The 4660 bump may be contami- 
nated by additional nebular emission from He II, 
He I, and forbidden Fe and Ar lines observation 
at medium or low resolution, the distinction be- 
tween the individual components of the bump is 
often difficult. This makes comparisons to the 



As Fig. All shows, the equivalent width of 
the 4660 bump is usually expected to be larger 
than that of 5808. In WC rich phases, however, 
VF(5808) can reach values comparable to that 
of 4660. Note that the appearance of the 4660 
bump before the WR phase (0.5 ^ t ^ 2.5 Myr) 
is due to stellar He II A4686 emission from Olf 
stars (cf. below). Compared to the recent mod- 
els of Cervino & Mas-Hesse (1994) ^aur predicted 
equivalent widths for the 4660 bump are larger 
by up to a factor of ~ 3 for low Z. The larger 
values are essentially due to the high mass loss 
evolutionary tracks adopted in the present work, 
which lead to larger WR/0 star ratios in bet- 
ter agreement with observations (see Maeder & 
Meynet 1994). 

For comparison we also plot the equivalent 
widths obtained using the bump luminosities of 
Smith (1991) given in Table ^ (dotted and short- 
dashed line in Fig. All). In this case larger EWs 
for the two WR bumps are expected as can be 
seen from a comparison of Tables | and | with 
p. Indeed the 4660 bump luminosity from Smith 
(1991), derived from He II A4686 of 5 LMC WN6- 
7 stars, represents an upper limit for all LMC 
WN stars (cf. Vacca 1992), even if the N iii/v A 
4640 contribution is included. It must, however, 
also be reminded that N lli/v A 4640 seems to 
be metallicity dependent (cf. § p. 2.1 ). We feel 
that, for most cases, our procedure of allowing 
for the strong variations of WC subtypes and 
including more recent calibrations for WN stars 
should yield more reliable predictions for the WR 
bumps than calculations based on the data of 
Smith (1991). 

Figure |A12| shows the relative WRbump/H/3 
line intensities of the two bumps. Given the 



® Their adopted WR bump luminosities are 1/4660 = 2 x 10'^^ 
erg for WN stars and L4650 = 6 x lO^*' erg s~^ for WCs, 
which differ by ~ 20-15 % from our values for WNL and 
WC4 stars. 
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steepness of the spectrum the maximum inten- 
sity ratio 5808/H/3 is considerably smaller than 
4660/H/3 in constrast to the maximum of the 
equivalent widths (cf. above). The strong in- 
crease of the peak 4660/H/9 with metallicity is in 
large part due to the much more rapid decrease 
of H/3 for large metallicities (cf. Arnault, Kunth 
k Schild 1989). 

6.2.2. WR humps at medium and high spectral 
resolution 



In Figure A13 we plot the time evolution of 
the equivalent widths of all individual WR lines 
(except Haand H/3; cf. § 6.2.4 ). The inclusion 



of 4686 A emission from Of stars, although its 
strength is considerably uncertain, leads to typi- 
cal values of VF(4686) ~ 1-2 A before the appear- 
ance of the first WR stars (t <^ 2-3 Myr). During 
the WN dominated phases the strongest lines are 
He II A 4686, N iii/v A 4640 and C iv A 5808 in 
decreasing order. The weakness of N iii/v A4640 
at subsolar metallicity is due to the metallicity- 
dependence of this line (SSM96). For solar and 
higher metallicities N iii/v A4640 is expected 
to have a strength comparable to He II A4686 
(Smith 1991, SMM96). In the WC-rich phase 
C lli/iv A 4650 is expected to dominate the 4660 
bump (I^(4650)/(T^(4640)-FVF(4686)) > 1). In- 
terestingly, for a given metallicity at Z > 0.001, 
the maximum equivalent width attained for the 
4650 bump. He II A 4686 and C IV A5808 are 
all of the same order. Since late WC stars are 
only expected at high metallicities (cf. § |3.2| ) the 
maximum equivalent width of C III A 5696 only 
reaches values above ~ 1 A for Z > 0.02. At sub- 
solar metallicities the equivalent width of He II 
A5412 is found to follow closely that of N iii/v 
A4640. In all cases VF(5412) is less than 2 A. 

Figures ^ (for Z=0.008, 0.004, and 0.001) 
and |AT| (for Z=0.02 and 0.04) show the ratios 
of the WR lines with respect to H/? as a func- 
tion of the H/3 equivalent width. These Figures 
involve purely observational quantities and can 
be used for direct comparisons with integrated 



spectra of a young burst population (with the 
aforementioned cautions in mind). Note that for 
all metallicities He II A4686/H/3 represents the 
maximum WR/HP intensity ratio obtained in 
a burst, although the maximum 4650 and 5808 
equivalent widths are comparable or even larger 
than 1^(4686). This is essentially due to stronger 
H/3 emission during the WC-rich phase at earlier 
ages, while 4686 emission from WN stars persists 
until later times where H/? has already consider- 
ably decreased. 

6.2.3. Nebular and stellar He II X4686 

The case of He II A4686 deserves closer exam- 
ination. As discussed in § || the models predict a 
non-negligible nebular He II A4686 emission dur- 
ing part of the WR-rich phase. It is therefore ex- 
pected that the 4686 A feature consists of both 
broad stellar and narrow nebular emission. Al- 
though observationally the distinction between 
narrow and broad components is not necessar- 
ily straightforward and has not generally been 
made, it is crucial for our understanding of exci- 
tation conditions and WR populations especially 
in low metallicity objects (cf. Garnett et al. 1991, 
and references therein; Schaerer 1996c, Izotov et 
al. 1996). A detailed discussion of the origin of 
nebular He II and comparisons with recent ob- 
servations will be presented elsewhere. Here we 
shall present only the main theoretical predic- 
tions. 



Figure A16 shows the intensity ratios A/H/3 of 
pure nebular He II A4686 (lower left) and both 
nebular and total WR+nebular 4686 emission 
(upper left) for metallicities 1/20 < Z/Zq < 1/5 
as a function of W(RP). For lower H/? equivalent 
widths (i.e. larger ages) 4686 is always dominated 
by stellar emission, and is therefore not shown 
here. The same also holds for solar and higher 
metallicities, as discussed by Schaerer (1996c). 
In the early phase of the burst nebular He II 
4686/H/3 reaches values up to ~ 1%, the max- 
imum value being nearly independent of metal- 
licity. Since nebular 4686 emission is related to 
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WC /wo stars and their hot WN star precursors 
(cf. § ^) whose population diminishes with Z, the 
duration of the phase with strong nebular emis- 
sion is also predicted to decrease. At Z =0.001, 
the current models predict that He II A4686 is 
dominated by stellar emission. 

In Fig. |A16| we have also plotted the predicted 
C III /iv A 4650 emission which is due to WC /WO 
stars. As mentioned ealier (§ p.2.1| ) this blend 
includes C III A4650, C IV AA4650,4658 and a 
small fraction of He II A 4686. Interestingly, at 
Z = 0.008 and 0.004, the predicted 4650 emis- 
sion has a strength similar to that of the total 
4686 emission. It should be remembered, how- 
ever, that at low Z, WC4 or WO stars are the 
dominant WC subtypes (Fig. [A^ ). In view of 
the strong variation of the adopted C IV emis- 
sion lines between WC4 and WO stars (see Table 
^ and present uncertainties in the assignment of 
WC/WO subtypes (cf. Smith k Maeder 1991, 
Kingsburgh, Barlow & Storey 1995), the predic- 
tions for 4650 remain somewhat uncertain. 

If nebular He II emission is indeed related to 
the WC/WO phase, one expects to detect other 
signatures for these stars. A nearly unambigous 
signature from WC/WO stars is C IV A5808 
whose strength is shown in Fig. A15. However, 



the same word of caution regarding the strength 
of 4650 also applies to the strength of 5808. The 
detection of O IV A3400, O VI A3811,34, or O V 
A5590 would be a clear signature of WO stars (see 
Kingsburgh, Barlow & Storey 1995). Test calcu- 
lations with the current models, however, predict 
very small equivalent widths for these features. 

6.2.4- Broad Ha and H/3 components 

The predicted emission in the He II A4861 + 
H/3 blend (denoted as 4861 in Tables || and|2[) and 
in the He II A6562 + Ha blend (denoted as 6560) 
is illustrated in Fig. A18| for all metallicities. For 
subsolar metallicities equivalent widths of ~ 5-25 
A for 6560 and 1-3 A for 4861 are predicted. For 
solar and larger metallicities stronger emission is 
expected. 



The line strengths predicted in this work do 
not include stellar absorption features. Such pre- 
dictions for young bursts have been given by 
Olofsson (1995a) who synthesized the major H 
and He absorption lines for instantaneous bursts 
at metallicities Z < Zq. For the ages shown in 
Fig. [A18 his predicted equivalent widths (absorp- 
tion) are: 2-4 A for Ha and 2-5 A for H/3 for the 
range of IMF considered. The strength of the 
absorption lines is fairly independent of metal- 
licity (Olofsson 1995a). (Note, however, that for 
unknown reasons the EWs of the Balmer lines 
predicted by Olofsson are systematically smaller 
than those of Diaz 1988.) The comparison with 
the results of Olofsson (1995a) shows that dur- 
ing the WR phase the equivalent width of the Ha 
emission from WR stars is larger than the equiva- 
lent width of the underlying absorption. The net 
result of the superposition depends, of course, on 
the relative widths of the two components. 

The strength of the stellar 6560 A emission is 
compared to that of the nebular emission in Fig- 
ure Al8| . For low metallicities the WR emission 
may reach several percent of the nebular flux, 
whereas much larger contributions are expected 
at solar and higher metallicity. This effect is bol- 
stered by both the increase of the WR population 
and the lower ionizing flux. As can be seen from 
Table ^ and | the relative 6560/4861 hue emis- 
sion in WR stars varies between 1.6 and 6. For 
an integrated population a typical value of this 
ratio is found to be ~ 4-5, which is larger than 
the value of Ha/H/3=2.86 for Case B recombi- 
nation. Thus, the relative contributions of WR 
emission to the nebular Ha and H/? lines may, 
in principle, affect reddening determinations de- 
rived from these lines. However, the uncertain- 
ties on the WR line emission in these lines are 
fairly large (cf. Tables ||, ^). 

Obviously the broad Ha and H/J emission 
components discussed here occur only during the 
WR phase, which is best detected by the presence 
of WR bumps that are uncontaminated (or only 
weakly contaminated) by nebular emission. In 
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this case the flux in the broad 4861 A emission 
provided by WNL stars is typicahy ~ 25 % of 
the flux in the total 4650 bump, while it is neg- 
ligible for WC stars (cf. Tables |, |). Additional 
broad emission which may be present (e.g. Roy et 
al. 1992, Izotov et al. 1996) cannot be attributed 
directly to WR stars. 

7. Models including massive close binary 
stars 

The main results derived from models in which 
binary stars are included according to the pre- 
scriptions in are illustrated in Figures |A19|and 
A20i for metalhcities Z=0.008 and 0.001. (The 
results for the remaining metalhcities show the 
same qualitative behaviour and therefore are not 
be reproduced here.) The top panel shows the 
relative fraction of WR stars of different sub- 
types. The use of WR/0 number ratios (as 
e.g. shown in Fig. |A2[) is avoided for two rea- 
sons. First, in the phase where the WR stars 
are formed through the binary channel (t ^ 5 
Myr) the number of single O stars is usually ei- 
ther low or possibly zero. Second, in order to 
compute such ratios one also needs to account 
for O stars formed from the mass-gainer in the 
binary systems (whose fate is not followed in the 
current calculations). In any case, as discussed 
earlier, it is advisable to directly use diagrams 
involving only observable quantities (cf. bottom 
panel). The middle panel shows the time evolu- 
tion of the equivalent widths of the strongest WR 
lines (He II A4686 and C IV A5808) as a function 
of time. The intensities of the same lines with 
respect to H/3 are plotted in the bottom panel as 
a function of the H/3 equivalent width. 

Accounting for the WR formation in massive 
close binary systems obviously extends the pe- 
riod where WR stars are present in instantaneous 
bursts (see top panels), since their progenitors 
can have masses smaller than the WR mass limit 
for single stars. For the binary models adopted 
here, WR stars can be present up to ages of ~ 
11-13 Myr. However, this period can be even 



longer (ages up to ^ 15-20 Myr) if the possible 
evolution of the secondary to the WR phase is 
also included in the model (Cerviho, Mas-Hesse 
& Kunth 1996, Vanbeveren, Van Bever & Bonder 
1997). In sharp contrast to burst models which 
include only single stars, the total duration of 
the WR phase is now essentially independent of 
metallicity. 

The time evolution of the relative fraction of 
the different WR subtypes (top panels) is eas- 
ily understood by combining the results given in 
Fig. A2 and Fig. 3 from de Loore &: Vanbeveren 
(1994), which give lifetimes in the WNL, WNE 
and WC phases function of the initial pri- 
mary mass. The formation of WR stars through 
the binary channel extends the flnal WNL phase 
of single stars (cf. §^) to later ages, leads to a 
second WC-rich period, and results in a signifl- 
cant phase dominated by WNE stars after about 
7-8 Myr. However, the following caveat regard- 
ing the WR subtype distribution in close binary 
systems should be kept in mind: As shown by 
Vanbeveren (1995) the predicted WC/WN num- 
ber ratio for low metallicity regions with constant 
star formation exceeds the observed values. One 
possible explanation is that the duration of the 
WC phase predicted from binary models is too 
long. Indeed the different observed WR fractions 
are well reproduced if RLOF leads to WR stars of 
mostly WNE type (see Maeder & Meynet 1994). 
Therefore, the relative distribution of WC and 
WNE subtypes in the binary phase should be 
considered fairly uncertain. 

The temporal evolution of the predicted equiv- 
alent widths of 4686 and 5808 (see middle panel) 
shows two main characteristics: 1) The phases 
where WR stars are formed predominantly through 
the single and the binary channels are well sep- 
arated. This behaviour, which is obviously true 
only for short bursts, is a direct consequence of 
the different progenitor masses from single and 
binary stars. 2) The expected maximum equiv- 
alent width of He II A4686 (C IV A5808) dur- 
ing the binary dominated phase is ~ 0.5-2 (2-4) 
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A, and is fairly independent of Z . 3) At low Z 
{Z ^0.001) the maximum equivalent widths of 
the WR features during the "binary dominated" 
phase become comparable or even stronger than 
those during the preceding "single star" phase. 

The bottom panel of Figs. [A19| and |A20| shows 
the predicted line ratios of 4686/H/3 and 5808/H/3 
as a function of the H/? equivalent width. As 
above the two phases where WR originate mostly 
from pure stellar wind mass loss ("single star 
channel") or from RLOF ("binary channel") are 
weU separated. For Z > 0.004 (Z=0.001) the 
binary-dominated phase corresponds to (H/3) ^ 
40 (100) A. The H/3 equivalent width (or other 
age indicator) should thus roughly provide a 
means to identify burst regions where WR forma- 
tion through RLOF may be of importance (see 
also Schaerer Sz Vacca 1996). Due to the strongly 
decreasing H/3 luminosity the WR/H/3 line ratios 
during the late binary phase can easily exceed 
the values during the early single-star dominated 
phase. The maximum line ratios of 4686/H/3 and 
5808/H/3 obtained for the lowest metallicities are 
typically 5 to 10 % of H/3 for the adopted binary 
fraction. These values might be overestimated 
by a factor of ~ 2 if at low Z the binary fre- 
quency is actually lower than the adopted value 

(cf. gmi) . 

Regarding the predictions including massive 
close binary systems we wish to stress again the 
numerous uncertainties mentioned in S; 2. 1.2. It 



is certainly true that our current knowledge of 
the evolution of binary stars is less accurate than 
that of single stars. In addition, in view of the 
simplified treatment adopted here for binary sys- 
tems we regard the detailed predictions as explo- 
rative and warn the reader from an overinter- 
pretation of this data. The general conclusions 
should, however, be fairly reliable. 



8. On determining the WR and O star 
content in H II regions and starbursts 

The absolute number of WR and O stars, 
as well as the relative WR/0 star number ra- 
tios in starburst regions, have been estimated 
with various techniques (e.g. Kunth & Sargent, 
1981; Arnault, Kunth & Schild 1986; Leitherer, 
1990; Mas-Hesse k Kunth 1991a, VC92, Vacca 
1994; Cerviiio & Mas-Hesse 1994; Meynet 1995; 
Schaerer 1996c). Knowledge of the absolute 
numbers of stars constituting a population can 
be used to determine quantities such as the star 
formation efficiency or rate in a region; relative 
WR/0 star number ratios can be used in con- 
junction with stellar evolution models to con- 
strain parameters such as the burst duration and 
the IMF slope, and to understand the variation 
in population parameters with metallicity. 

Rough estimates of the WR and O star con- 
tent in a starburst region can be obtained in a 
straightforward way. However, as demonstrated 
by Schaerer (1996c), more accurate techniques 
are required to obtain reliable constraints on 
both the O star content (cf. Vacca 1994) and the 
IMF (cf. Meynet 1995; Contini, Davoust k Con- 
sidere 1995). In the following we will therefore 
review the various methods employed and discuss 
procedures that should yield the most accurate 
results. 

8.1. Determination of the number of O 
stars 

Usually the number of O stars present in an 
integrated population is estimated from compar- 
ing 1) the total UV continuum luminosity with 
the predictions of evolutionary synthesis mod- 
els (e.g. Leitherer k Heckman 1995), or 2) the 
observed strength of the H recombination lines 
with predictions of the ionizing fluxes from stel- 
lar atmosphere models (e.g. Osterbrock 1989). 
Here we shall concentrate on the latter proce- 
dure, which has been discussed extensively by 
Vacca (1994, and references therein). 
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8.1.1. Using H recombination lines 

This method is based on the premise that the 
number of ZAMS OV stars can be directly re- 
lated to the observed H/? luminosity (cf. Vacca 
1994). The ionizing flux (in terms of the total 
number of H-ionizing photons, Qg"*'^^) required 
to produce the H/3 line luminosity L(H/3) is given 

by 

.Total _ aB(H°)AH/3 



-L{m 



(7) 



(cf. Eq. Q with = 1). Note that Eq. ^ assumes 
a dust-free and ionization-bounded H II region 
(Case B recombination). If the stellar popula- 
tion is the only source of ionizing photons, Qo"*^' 
can be expressed in terms of the number of so- 
called "equivalent" stars of a given subtype (here 
we use Nqjy to denote the number of "equiva- 
lent 07V" stars) responsible for producing the 
ionizing luminosity. 



Q 
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07V 



(8) 



where Qq"^^ is the Lyman continuum luminosity 
of an individual 07V star. As shown by Vacca 
(1994) the number of ZAMS O V stars is then 
related to iVo7V 



Kyy/NoY, 



(9) 



The number of ZAMS OV stars can thus be eas- 
ily derived from the observed H/3 luminosity us- 
ing the extensive tables of Vacca (1994) which 
provide ijq for different assumptions of the IMF 
slope, mass limits and metallicity. For a Salpeter 
IMF and an upper mass cut-off of 120 Mq e.g., 
rjo is found to be close to unity. 

However, for comparisons with evolutionary 
models and other applications it is crucial to dis- 
tinguish further between ZAMS OV stars and O 
stars in general. Indeed, for instantaneous bursts 
the number of ZAMS OV stars obtained from 
the above procedure approximates the total num- 
ber of O stars (of all luminosity classes) only for 
ages less than ~ 2-3 Myr, as demonstrated by 
Schaerer (1996c). For larger ages, A^'ov derived 
with the method of Vacca (1994) generally over- 
estimates the true number of O stars, if one does 
not appropriately reduce the upper mass limit. 
Nevertheless, the procedure of Vacca (1994) can 
be easily generalized to account for the tempo- 
ral evolution of the total ionizing luminosity and 
thus to derive the total numbers of O stars. Ac- 
cordingly, Schaerer (1996c) introduced the time 
dependent value of r/o: 



^0 (t) 



j;^;z<^{M)Qo {M,t)dM 



(12) 



where 



i07V r^"p 



(10) 



represents the IMF averaged ionizing Lyman con- 
tinuum luminosity of a ZAMS population nor- 
malized to the output of one "equivalent" 07V 
star (see Vacca 1994). In Eq. 10, <1> represents 



the IMF with the mass limits Miow and M^p , and 
Qq^^^{M) is the Lyman continuum luminosity 
as a function of stellar mass for ZAMS stars. The 
number of ZAMS OV stars is thus given by 



iVov 



Total 



(11) 



where Qq {M, t) is the ionizing luminosity of a 
star of initial mass M at time t. Note that for 
the integral in the denominator, the identical def- 
inition for O stars is used as in the evolutionary 
synthesis models. With Eq. |l^ the total number 
of O stars is now given by 



No 



07V 



(13) 



The temporal evolution of ??o(^) for our "stan- 
dard" instantaneous burst models is shown in 
Fig. A21 for all metallicities considered in this 
work. As expected, the ZAMS values {r]o{t = 0)) 
agree well with the results of Vacca (1994) using 



the same IMF. Fig. A21 shows that r/o remains 
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fairly constant within the first ~ 2 Myr, and 
then decreases due to the evolution of the O star 
population and the subsequent disappearance of 
the most massive stars (cf. Schaerer 1996c). For 
subsolar metallicities a short-lived peak of 779 
is found at the beginning of the brief WR-rich 
phase (t ~ 3 Myr) due to the large average Ly- 
man continuum luminosity of WR stars at this 
time (cf. Fig. |A6| ). The second, brief increase in 
?7o seen in several of the panels occurs when the 
number of O stars (as defined by the lower limit 
on Tgff ) drops to zero, and hence the denominator 
in Eq. ^ approaches zero. Obviously, the defini- 
tion of rjo becomes meaningless at this time; the 
ionizing flux of the population will be provided 
by stars of B type or later and WR stars which 
may still be present in some cases {Z > 0.02). 

Although more general than the formulation 
of Vacca (1994), the explicit time dependence of 
770 (t) introduces two new variables that must be 
determined before this quantity can be used to 
determine the massive star population in a re- 
gion. In addition to the assumptions about the 
IMF (slope and mainly upper mass limit) and 
metallicity, the value of r]o{t) also depends on 
the age of the stellar population and the star 
formation history, for which we have only con- 
sidered the case of an instantaneous burst here. 
In most cases, the observed H/3 equivalent width 
can be used to estimate the population age (see 
Fig. ^^) . The star formation history is somewhat 
more problematic, as this is usually one of the pa- 
rameters trying to be determined. Note that the 
procedures and limitations discussed here also 
apply in principle to methods which estimate the 
O star content from the UV continuum luminos- 
ity. However, these methods tend to be less sen- 
sitive to variations in the population parameters 
and avoid other potential difficulties affecting the 
nebular lines (photon leakage etc.). 

8.1.2. Using He recombination lines 

An equivalent of the above procedure can be 
formulated using the He I A4471 line, which relies 



on Qi the total He I ionizing luminosity emitted 
by the stellar population(cf. Vacca 1994). How- 
ever, since stars with Tgfj ;^ 40000 K provide 
enough photons above 24.6 eV to fully ionize 
helium, the ratio of the ionic fractions of He I 
and H I in the nebula becomes independent of 
Teff (see e.g., Doherty et al. 1995). In this case 
He I does not provide any additional information 
compared to H recombination lines. In partic- 
ular He I A4471/H/3 will be insensitive to vari- 
ations in the IMF (both slope and upper mass 
limit) as long as enough hot stars are present in 
the population to keep the "equivalent" effective 
temperature above ~ 40000 K. 

The situation is even more complicated for 
other He I lines, e.g. the He I A2.06 /um line, 
whose predicted strength is dependent on as- 
sumptions about the structure of the H II re- 
gion (cf. Shields 1993). Lines resulting from re- 
combination of He"^"^ are rarely detected in H II 
regions and their origin is not fully understood 
(cf. Schaerer 1996c and references therein). Dur- 
ing the pure O star phase of a burst (cf. § ^) the 
He II A4686/H/3 is most likely too weak to be de- 
tected. In later phases WR stars (or other poorly 
understood processes) dominate the formation of 
nebular He II lines. See § ^ for more details. 

8.2. Determination of the number of WR 
stars 

The total number of WR stars, and their rel- 
ative subtype distribution, can be easily deter- 
mined from the observed luminosity in a given 
broad WR line in the integrated spectrum of a 
starburst region (e.g., Osterbrock &: Cohen 1982; 
Kunth & Sargent 1983; VC92). The rehability of 
this method obviously rests on the constancy of 
the line luminosity for all WR stars (cf. § 2.2.1) 
and a correct measurement of the desired WR 
component of the line, which is often a compos- 
ite of different WR lines (see § p. 2.1 ) or blended 
with nebular emission (e.g. [Fe III] A4658, [Ar IV] 
AA4711, 4740). 
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8.3. Determination of WR/O number ra- 
tios 

Two different methods have been developed 
to derive WR/O star number ratios from op- 
tical spectra of integrated populations. While 
the method of Arnault, Kunth & Schild (1986) 
uses the flux in the entire WR bump, the more 
recent work of VC92 relies only on the flux in 
the broad He II A4686 line. In the following we 
will briefly discuss and compare these two meth- 
ods, point out their respective advantages, and 
present some improvements on these techniques. 

8.3.1. From He II \4686 and H/3 

Following VC92, we can write the ratio of the 
number of WR stars to the number of reference 
07V stars as 



Total 



(14) 



where is the average Lyman continuum lu- 
minosity per WR star. The total number of WR 
stars is given by the luminosity of the broad WR 
line, 

AwR = -^'obs(AwR)/-^^WR(AwR)> (15) 

where AwR represents one particular WR line 
(e.g. He II A4686), and -^^wr('^wr) is the line 
luminosity for a single WR (e.g. WN) star. In- 
serting Eqs. 0, ^ and 15 in Eq. ^ one finds. 
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(16) 



For the nebular conditions assumed here (cf. § 
P) aB(H0)AH/3/(aff^M = 2.1 x lO^^ erg-^. In 
addition to the assumptions stated earlier re- 
garding the presence of dust and leakage of ion- 
izing photons (cf § 8J), Equation |l^ assumes 



that the ionization is provided completely by the 
population of Nq O stars and A^wr WR stars 



with a Lyman continuum luminosity Qq and 
Qo^^, respectively. Furthermore it is assumed 
that only stars of the generic subtype "WR" con- 
tribute to the WR line in consideration. Vacca 
& Conti (1992) considered the case where only 
WNL stars are present. They adopted the values 
of Lwnl(A4686) = 1.7 x 10^6 erg s'^, Q^^^ = 
1.70 X 10^9 s-\ and Q^'^^ = 1.01 x 10^^ s-\ in 
Eq. 15 of VC92 for %(*) = 1- 

We will now try to estimate the accuracy of 
the procedure of VC92, or more generally of 
Eq. |l^. In addition to the age of the popu- 
lation and the star formation history (needed 
for r/o; § |8.1D , use of Eq. ^ requires an esti- 
mate of Qo^^/ Qq^^ , the ratio of the average Ly- 
man continuum luminosity emitted per WR star 
to the Lyman continuum output of the equiva- 
lent 07V star. As shown in Fig. large vari- 
ations of Qq^^ are expected during the evolu- 
tion of a burst. Values of logQo^^ > 49.05 
(i.e. Q^^/Q§^ > 1.68) yield larger WR/O ra- 
tios than those derived using Eq. 15 of VC92. 
However, for the typical values of He II A4686/H/?< 
0.1 observed in starburst regions, an average 
value of below the value used by VC92 

yields negligible changes in Awr/A^o- 

Equation ^ accounts for the presence of one 
dominant WR subtype, e.g. WNL stars as as- 
sumed by VC92. However, as discussed in § 
^, variations in the subtype distribution are ex- 
pected clS 3j burst population ages. Although an 
appropriate WR line, whose flux is due primarily 
to WR stars of a certain subtype (i.e. 4686 A by 
WNL stars, or 5808 A by WC stars) could easily 
be chosen in the case of a mixed population, for 



Eq. |16|to be valid the representative WR subtype 
must also contribute the bulk of the ionizing flux 
from the total WR population. 

In conclusion we stress again (cf. Schaerer 
1996c) that variations of ??o(t) and Qq^^ compli- 
cate the derivation of accurate WR/O number 
ratios from He II A4686/H/3 using the method 
devised by VC92. For detailed comparisons with 
evolutionary models and reliable estimates of 
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other parameters like the IMF slope etc., it is 
preferable to use models such as those developed 
in the present work, which allow a direct compar- 
ison of all relevant observational quantities (see 
Sects. and thereby circumvent the difficulties 
discussed above. 

8.3.2. From the "WR-bump" and H/3 

From the first synthesis models of WR and O 
star populations Arnault, Kunth & Schild (1986) 
proposed an approximate fit relation to derive an 
estimate of the total WR/0 ratio (for all WR 
subtypes) from the ratio of the WR bump lumi- 
nosity to H/3. To the level of accuracy consid- 
ered in their work the results for a given IMF 
slope were found to be fairly insensitive to the 
star formation scenario and metallicity. Despite 
this insensitivity, to the best of our knowledge the 
fit relation proposed by Arnault, Kunth & Schild 
(1986) has not been used in the literature to es- 
timate WR/0 ratios. With the new models pre- 
sented in this work, which rely on more complete 
input physics, we have reexamined the results of 
Arnault, Kunth & Schild (1986). To calculate the 
total emission in the WR bump, we sum the in- 
dividual contributions of N iii/v A4640, C iii/iv 
A4650 and He II A 4686 from WN and WC stars. 



Figure |A22| shows the WR/(WR-F0) ratio as a 
function of the predicted WR-bump/H/3 ratio for 
all instantaneous burst models (thin lines). 

As anticipated from the results of Arnault, 
Kunth & Schild, the data show reasonably well- 
defined relation extending over at least 2 orders 
of magnitude in WR-bump/H/?. Our data can 
be fitted by 



log 



WR 



(WR + O) 



(-0.11±0.02)-h 



(0.85 ± 0.02) log ( WRbump^ (-^7) 

(thick solid line labeled SV) with a rms of 0.22 
dex. As discussed below (cf. § 5.2.1 ), the de- 
tailed predictions show notable differences if we 
adopt the WR bump luminosities of Smith (1991, 



cf. Table ^). The average fit relation is, how- 
ever, only slightly changed by this assumption. 
For comparisons the thick dotted line (labeled 
AKS86) shows the fit relation from Arnault, 
Kunth & Schild (1986). Although our detailed 
predictions reveal large differences with the re- 
sults of Arnault, Kunth & Schild (1986), the de- 
rived fit relations agree relatively well. 

We have also compared the model predictions 
with the WR/(WR-|-0) ratio derived from Eq. |l^ 
assuming %(*) = 1, Qq^^/Qq'^^ = 1-68 and two 
values of the WR-bump luminosity L-wr(AWR) = 
1.7 X 10^^ erg s~^ (thick long-dashed line) and 
5. X 10'^^ erg s"^ (thick short-dashed line) respec- 
tively. The long-dashed line thus corresponds 
to Eq. 15 of VC92, while the short-dashed line 
stands for the upper limit of the WR-bump lu- 
minosity (WC stars from Smith 1991). The com- 
parison of both curves with the model predictions 
confirm the finding that in general the method 
of VC92 overestimates the WR/0 ratio. The 
largest differences are obtained for large values 
of WR-bump/II/3, which are expected at late 
epochs in bursts. In this case the difference is 
due to the assumption of an unevolved ZAMS 
population, i.e. rjQ = const (cf. above). 

We conclude that Eq. [l^ is useful for estimates 
of the total WR/(WR-|-0) ratio from relatively 
low resolution integrated spectra of starburst re- 
gions. As a word of caution it must be mentioned 
that such a relation does not hold for the He II 
A4686/H/3 (or even for C IV A5808/H/?) as can 
be seen from Figs. and A14. The fit rela- 



tion also requires that all possible contamination 
from nebular lines (cf. § ^]^) to the WR bump be 
excluded before it can be used to determine the 
WR/0 star number ratio. In general it may not 
be possible to do this. 

8.3.3. From UV spectra 

Since the luminosity of WN and early WC 
stars in the He II A1640 line is quite similar (see 
Tables II), and since O stars are the dominant 
contributors to the UV continuum, the equiva- 
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lent width of He II A 1640 is expected to be closely 
related to the ratio of the total number of WR 



stars to O stars. Indeed, as Fig. [A23| shows, our 
models predict a relatively well-defined relation 
between 1^(1640) and WR/(WR+0), which can 
be fitted by 

WR 



log 



-1.31 ±0.03) + 



(WR + 0). 
(1.52 ± 0.05) log (EW(1640)) (18) 

with a rms of 0.29 dex (thick solid line in Fig. [A23D , 
where VF(1640) denotes the equivalent width of 
the broad stellar He II A1640 emission in A. Note 
that in this fit we have excluded the points with 
WR/(WR+0)=1, corresponding to bursts with 
ages t ^ 3-4 Myr at Z > 0.02. 

Equation |l^ can be used to obtain estimates 
of the total WR/0 ratio from a pure equiva- 
lent width measurement in cases where the rest- 
frame UV is the only accesible wavelength range 
(e.g. high redshift galaxies), or if the optical spec- 
trum is strongly contaminated by old popula- 
tions. As long as there is no emission compo- 
nent from an active galactic nucleus, the He II 
A1640 line also has the advantage of being inde- 
pendent of nebular lines, which may be subject 
to additional complications, as mentioned above. 
For more precise determination of the WR and 
O star content, comparisons with the detailed 
model predictions discussed in § |6| should be pre- 
ferred over the use of Eq. 18. 



9. Summary 

We have constructed new evolutionary synthe- 
sis models for young starbursts using the latest 
stellar evolution tracks, theoretical spectra gen- 
erated with the most complete input physics, and 
a compilation of observed emission line strengths 
for WR stars. To provide predictions for various 
stellar and nebular emission line strengths, our 
models account explicitly for the the contribu- 
tions of WR stars of different subtypes. We have 
also explored the consequences of WR formation 
through mass transfer in close binary systems. 



9.1. Predicted observational features 

The models presented here, and their predic- 
tions for various emission line strengths, can be 
used to derive massive star populations in star- 
burst regions. Comparisons between observed 
and predicted line strengths allow the derivation 
of several properties characteristic of a starburst 
(e.g. burst age and duration, IMF, mass). In par- 
ticular, we synthesize seven different WR emis- 
sion lines, including one UV line. We also predict 
the strength of the various nebular and stellar 
components contributing to the so-called WR- 
bump (A ~ 4650 A). These lines, whose emission 
strengths vary considerably among the various 
WR subtypes, provide a means of deriving the 
numbers of both WN and WC stars in a star- 
burst region. The model predictions for the neb- 
ular lines provide a measure of the O star popu- 
lation. In addition, our models provide the first 
quantitative predictions regarding a possible ori- 
gin for nebular He II emission which is frequently 
observed in low metallicity H II regions (e.g. Gar- 
nett et al. 1991). We also estimate the contribu- 
tion of WR emission to broad components seen 
underlying Ha and H/3 lines in the spectra of H 
II regions (e.g. Roy et al. 1992). Taken together 
these diagnostic features should allow the most 
accurate determination of massive star popula- 
tions in young bursts to date. 

9.2. Stellar evolution in different envi- 
ronments and tests of the models 

Models of massive star evolution and atmo- 
spheres have not been extensively tested in re- 
gions outside of the solar vicinity. It is particu- 
larly important to understand the role metallicity 
plays in determining the characteristics of mas- 
sive stellar atmospheres and evolution. Our mod- 
els can be used to study the evolution of massive 
stars in environments which are not accessible 
"locally" (e.g. in very low or very high metal- 
licity regions), and therefore they can be used 
to test the models of stellar evolution and stel- 
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lar atmospheres under diverse and extreme con- 
ditions. For example, by comparing our model 
predictions with the parameters derived from in- 
tegrated spectra of starforming clusters whose 
stellar content and burst properties can be de- 
termined by other means, the effects of metallic- 
ity and environment on massive star properties 
and populations can be probed. Although such 
studies have only recently begun for LMC and 
SMC clusters whose content is known from stel- 
lar censuses (i.e., from imaging and direct count- 
ing; Vacca et al. 1995; Vacca &: Conti 1997), they 
will serve as touchstones for evolutionary synthe- 
sis models. Similar studies can be carried out for 
other star-forming regions where individual stars 
can be resolved with HST. It is particularly inter- 
esting to extend such studies to extremely metal 
poor galaxies (e.g., I Zw 18, Hunter & Thronson 
1995). 

There are several observational studies that 
can be carried out to test the validity of our mod- 
els: 1) Determination of the frequency of WR 
rich starbursts and the frequency of WC detec- 
tions among WR galaxies, as first suggested by 
Meynet (1995) (see §^). 2) Examination of the 
predicted relation between nebular He II A4686 
and a population of WC/WO and/or hot WN 
stars (cf. 3) Searches for WR signatures 

in old bursts, which may be indicative of WR 
stars formed through mass transfer in close bi- 
nary systems (§0). The first investigations along 
these lines have recently begun (Schaerer 1996c, 
Schaerer et al. 1997). 

Of particular interest in this regard is the anal- 
ysis of the optical spectra of the recent samples 
of low metallicity extragalactic H II regions, used 
for the determination of the primordial He abun- 
dance (e.g. Izotov, Thuan & Lipovetsky 1997; 
Olive, Skillman &: Steigman 1997). These objects 
are ideal for testing the evolutionary and atmo- 
spheric models in low Z environments. This is 
nicely illustrated by the recent detection of WC 
stars in I Zw 18 (Izotov et al. 1997a, Legrand et 
al. 1997). Studies of stellar populations in Seyfert 



galaxies and AGN (cf. Heckman et al. 1997) may 
also provide new insights into stellar evolution in 
metal rich environments. 

9.3. Future applications 

With the results from such analyses of stel- 
lar populations in hand, additional studies of the 
starburst phenomenon can proceed quickly. For 
example, based on the modelling of the WR and 
O populations in large numbers of starburst re- 
gions, possible variations from region to region in 
the the slope of the initial mass function for mas- 
sive stars can be investigated (cf. Schaerer 1996c; 
Mas-Hesse & Kunth 1997). Also the detection- 
frequency of WC stars in WR galaxies could pro- 
vide limits on the burst duration (cf. §^). Fur- 
thermore, the ejecta from a large population of 
WR stars may have a significant impact on the 
chemical evolution of the host galaxy (cf. Pagel 
et al. 1992; Esteban & Peimbert 1995; Kobul- 
nicky et al. 1997). With detailed knowledge of 
the WR and O star populations in starburst re- 
gions, the possible links between the massive star 
content and spatial variations in metallicities can 
be investigated. 

We envision that the present models will be 
used primarily to interpret the spectra of WR 
galaxies. These objects provide a window on the 
processes occurring during the youngest phases 
of a starburst (e.g. VC92; Contini, Davoust & 
Considere 1995; Conti 1996, Leitherer et al. 1996; 
Schaerer 1996c; Gonzalez-Delgado et al. 1997, 
Mas-Hesse & Kunth 1997; Schaerer et al. 1997). 
As such, a study of their properties has rele- 
vance to our understanding of objects ranging 
from nearby H ii regions to high redshift galax- 
ies. 

Future studies of starburst regions will move 
beyond the determination of the stellar popula- 
tions and will focus on broader aspects of the 
starburst phenomenon, such as the chronology 
of starburst events, the initial mass function and 
its possible dependence on environmental effects, 
mixing processes in the ISM and the chemical 
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evolution of young systems. Obtaining a com- 
plete and consistent picture of a starburst region, 
incorporating both the stellar content with its ra- 
diative, mechanical and chemical output and the 
evolution of the ISM of the host galaxy would be 
a major goal of these efforts. Our models should 
provide a first step toward this goal. 

We thank Paul Crowther, Laurent Drissen and 
Alex de Koter for kindly providing us with mea- 
surements from their spectra. We also thank 
Paul for answering a number of questions regard- 
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earlier version of the manuscript. DS is sup- 
ported by the Swiss National Foundation of Sci- 
entific Research and acknowlegdes partial sup- 
port from the Directors Discretionary Research 
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in the form of a fellowship from the Beatrice Wat- 
son Parrent Foundation. 

A. Appendix 

The results from the models presented in this 
work are available in electronic format. Addi- 
tional data such as spectral energy distributions 
and integrated colors are also provided. Further 
variations of the major input parameters (IMF 
slope etc.) are also considered. The models are 
found on the Web at 

jittp : / / www . stsci . edu/ ftp/ science/starburst 
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Fig. Al. — Luminosity in the He II 4686 line vs. 
the bolometric luminosity for WNL stars, from 
the measurements and analysis by Crowther &; 
Dessart (1997). Open circles denote LMC stars; 
filled squares are Galactic stars. Note the huge 
variation in line luminosity across the WNL sub- 
class. 
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Fig. A2. — Predicted Of/0 and WR/0 star num- 
ber ratios clS Si function of age for all metallici- 
ties. All models are calculated for an instanta- 
neous burst with a Salpeter IMF. Shown are the 
total WR/(WR+0) (solid), the WNL(WR+0) 
(dashed), the WC/(WR+0) (dashed-dotted) , 
and the Of/0 number ratios. Note the change of 
the ordinate between the top, middle, and lower 
panels. 
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Fig. A3. — Temporal evolution of the IMF- 
averaged WC subtype for all metallicities 
(Z=0.040: dashed-dotted, 0.020: short dashed, 
0.008: solid, 0.004: long dashed, 0.001: dotted). 
WO stars are designated by a WC3 subtype 



Fig. A4. — Temporal evolution of the ionizing 
photon flux of a stellar population (normalised 
to a total mass of 1 Mq) for all metallicities. 
Top panel: Lyman continuum flux, Middle: Flux 
in the He i continuum. Lower panel: flux in the 



He II continuum. Same symbols as in Fig. A3 
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Fig. A5. — Comparison of predicted ionizing 
fluxes using difl'erent atmosphere models. Shown 
is the standard model for Z=0.020. Dotted: us- 
ing plane parallel, LTE, line blanketed models 
of Kurucz (1992) only; Dashed: same as dotted 
but using spherically expanding non-LTE mod- 
els for WR stars (Schmutz, Leitherer & Gruen- 
wald, 1992); Solid: same as dashed but using 
spherically expanding, line blanketed, non-LTE 
CoStar models for O stars. The top, middle, and 
bottom panel show Qi,and Q2 respectively 
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Fig. A6. — Time evolution of the average Lyman 
continuum luminosity provided per WR star in 
instantaneous burst models. The thin solid hor- 
izontal line indicates the value of Q^^^ from 
Vacca & Conti (1992). Same symbols as in 



Fig. A3 
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Fig. A7. — Time evolution of the H/3 equivalent width for an instantaneous burst at the given metallicity. 
Same symbols as in Fig. |A3| . Left panel: Logarithm of M^(H/3) for ages from to 10 Myr. Right panel: 
Evolution of VF(H/3) during the WR rich phases of the burst. Horizontal lines on the margins show the 
maximum/minimum l/^^(H/3) during the WR phase 
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Fig. A8. — Predicted nebular emission line ratio 
of He II A4686/H/3 function of age. Same 



symbols as in Fig. A3. See Figs. A13 and A16 
for the total WR and nebular 4686 emission 
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Fig. A9.— Lower panels: total WR/(WR+0) 
(solid lines), WNL/(WR+0) (dashed), and 
WC/(WR+0) (dashed-dotted) number ratios 
as a function of time for Z=0.004 (left) and 
Z=0.001 (right). Upper panels: Nebular He II 
A4686/H/3 as a function of age for Z=0.004 (left) 
and Z=0.001 (right). The dotted hues show the 
predicted He II A4686 /H/5 when WC stars are ex- 
cluded. This shows that both WC and (hot) WN 
stars contribute to the nebular He II emission 
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Fig. AlO.— Evolution of the He II A1640 equiv- 



alent width. Same symbols as in Fig. AS 




Fig. All. — Evolution of WR bump equivalent 
widths for all metallicities. Solid lines: 4660 
bump including N iii/v A4640, C iii/iv A4650, 
and He II A4686. Long-dashed: C IV A5808. 
Both following the prescriptions in Tables || and 
0. Dotted and short-dashed lines show the WR 
bump equivalent widths using the prescriptions 
of Smith (1991; cf. Table pt) 
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Fig. A12. — Evolution of WR bump intensities 
relative to H/? for all metallicities. Solid: 4660 
bump, long-dashed: C IV A5808, computed as in 



Fig. All. The less accurate values using the pre- 
scription of Smith (1991) is not shown here. Cer- 
tain cases show non-uniqueness of line intensities 
when plotted against W(RP) (e.g. for Z=0.001, 
but cf. also other Figs.). This behaviour is due 
to the slight increase of the H/3 with time at the 
beginning of the WR phase (see Fig. [AT] ), which 
is in turn related to the contribution of WR stars 
to the Lyman continuum flux 
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Fig. A13. — Predicted equivalent widths of optical WR lines for all metallicities. Left panels: He IIA4686 
(solid), C iii/iv A4650 (dotted), N lli/v A4640 (short dashed-dotted) . Right panels: C IV A5808 (long 
dashed), C III A5696 (short dashed). He IIA5412 (long dashed-dotted) 
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Fig. A14. — Predicted intensity ratio /(A)//(H^) 
for WR lines at Z=0.008 (solid), 0.004 (long- 
dashed), and 0.001 (dotted) 
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Fig. A15.— Same as Fig. ^l|for Z=0.02 (short- 
dashed) and 0.04 (dashed-dotted) 
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Fig. A16.— 
Predicted intensity ratios I{X)/ I(H/s) for nebular 
and WR lines at A ~ 4650-4700 A for Z=0.008 
(solid), 0.004 (long-dashed), and 0.001 (dotted) 
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Fig. A17. — Predicted equivalent widths of Ha 
(upper panel) and H/3 (lower panel) blends for all 
metallicities. Symbols as in Fig. 



46 



0.040 
0,020 
0,008 
0,004 
0,00 




Fig. A18. — Ratio of broad WR emission to neb- 
ular Ha line for all metallicities. Symbols as 



in Fig. |A^. The relative contribution to H/3 is 
nearly identical to the Ha contribution (see § 
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Fig. A19. — Model predictions including massive 
close binary stars with a binary fraction / = 0.2 
for Z=0.008 (instantaneous burst). Top panel: 
Relative fraction of WR stars of different sub- 
types. Middle: Temporal evolution of predicted 
equivalent widths for He II A4686 (solid) and 
C IV A5808 (long-dashed). The arrow indicates 
the beginning of the phase where WR stars are 
predominantly formed through the binary chan- 
nel. Lower panel: Line ratios A/H/3 as a function 
W{Ilf3). Same symbols as middle panel. The 
vertical dotted lines illustrates the clear separa- 
tion of the "single star" and "binary" phases in 
this diagram 
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Fig. A20.— Same as A19 for Z=0.001 
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Fig. A21. — Evolution of r/o(t) (defined by Eq. 
^) in standard burst models for metallicities 
from Z=0.040 to 0.001 
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Fig. A22.-- Number ratio WR/(WR+0) as a 
function of the intensity ratio of the WR-bump 
over H/3 obtained from instantaneous burst mod- 
els at all metallicities (thin lines using same sym- 
bols as in Fig. The thick sohd line (SV) 
shows the fit to our models (cf. Eq. |l^). Thick 
dashed: relation from Arnault, Kunth & Schild. 
Thick long-dashed: relation from VC92 (their 
Eq. 15 assuming rj = 1.). Thick short-dashed: 
same as long-dashed but assuming a larger aver- 
age line luminosity of log(4650-4686)=36.7 
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Fig. A23.— Number ratio WR/(WR-F0) as a 
function of the He II A 1640 equivalent width 
(emission, in A) obtained from instantaneous 
burst models at all metallicities (thin lines us- 
ing same symbols as in Fig. [A3|). The thick solid 



line shows the fit to our models (cf. Eq. 18) 
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